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Abstract. We consider an incompressible viscous flow without surface tension in a finite- 
depth domain of three dimension, with free top boundary. This system is governed by a 
Naiver-Stokes equation in a moving domain and a transport equation for the top boundary. 
Traditionally, we consider this problem in Lagrangian coordinate and perturbed linear form. 
In [1], I. Tice and Y. Guo introduced a new framework using geometric structure in Eulerian 
coordinate to study both local and global wellposedness of this system. Following this path, we 
extend their result in local wellposedness from small data case to arbitrary data case. Other 
than the geometric energy estimate and time-dependent Galerkin method introduced in [1] , we 
■ utilize a few new techniques: (1) using parameterized Poisson integral to construct a nontrivial 

transform between fixed domain and moving domain; (2) using bootstrapping argument to 
prove a comparison result for steady Navier-Stokes equation for arbitrary data of free surface. 
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1. Introduction 

cl': 

1.1. Problem Presentation. We consider a viscous incompressible flow in a moving domain. 

jh- (1.1.1) n{t) = {yeExR\-b{y 1 ,y 2 )<y 3 < V {y 1 ,y 2 ,t)} 

Here we can take either E = R 2 or £ = (LiT) x {L 2 T) for which T denotes the 1-torus and 
L X ,L 2 > the periodicity lengths. The lower boundary b is fixed and given satisfying 

(1.1.2) < b(yx,y 2 ) <b = constant and b G C°°(£) 

When £ = R 2 , we further require that b{y\, y 2 ) — > bo = positive constant as \yi \ + \y 2 \ — > oo and 

b-b G # s (£) for any s > 0. When £ = (LiT)x(L 2 T), we just denote b = l/2(max{6(j/i,y 2 )}+ 

mm{b{y%,y 2 )}). It is easy to see this also implies b — bo G _fF(£) for any s > since £ is a 

bounded domain. We denote the initial domain Q(0) = Oo- For each t, the flow is described by 

| velocity and pressure (u,p) : Q(t) i— )■ i? 3 x R which satisfies the incompressible Navier-Stokes 

O ■ equation 
(N 



in 



X 



(1.1.3) 



dtu + u ■ Vit + Vp = fiAu 


in 


fi(t) 




V • u = 


in 


fi(t) 




{pi — (jM{u))v = gr\v 


on 


{Y3 = 


?y(yi,y2,t)} 


u = 


on 


{Y3 = 


-b(yi,y 2 )} 


dtV = u 3 - u\d yx r} - u 2 d y2 rj 


on 


{Y3 = 




u{t = 0) = Uq 


in 






rj{t = 0) = r/o 


on 


£ 





for ^ the outward-pointing unit normal vector on {2/3 = rj}, I the 3x3 identity matrix, (Bii)y = 
diiij +djUi the symmetric gradient of n, 5 the gravitational constant and /U > the viscosity. As 
described in [1], the fifth equation in (1.1.3) implies that the free surface is convected with the 
fluid. Note that in (1.1.3), we have make the shift of actual pressure p by constant atmosphere 
pressure p atm according to p = p + gy 3 - p atm . 

We will always assume the natural condition that there exists a positive number p such that 
r /o + ^^/°>Oon£, which means that the initial free surface is always strictly separated from 
the bottom. Also without loss of generality, we may assume that p = g = 1, which in fact will 
not infect our proof. In the following, we will use the term "infinite case" when £ = R 2 and 
"periodic case" when £ = {L\T) x {L 2 T). 
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1.2. Previous Results. Traditionally, this problem is studied in several different settings ac- 
cording to whether we consider the viscosity and surface tension and the different choices of 
domain. Under each setting, many different kinds of framework have been developed. 

In the viscous case without surface tension, the local wellposedness of equation (1.1.3) was 
proved by Solonnikov and Beale. Solonnikov [11] employed the framework of Holder spaces 
to study the problem in a bounded domain, all of whose boundary is free. Beale [2] utilized 
the framework of L 2 spaces to study the domain like ours. Both of them took the equation 
as a perturbation of the parabolic equation and make use of the regularity results for linear 
equations. Abel [20] extended the result to the framework of L p spaces. 

Many authors have also considered the effect of surface tension, which can help stabilize the 
problem and gain regularity. However, most of these results are related to global wellposedness 
for small data. 

Almost all of above results are built in the Lagrangian coordinate and ignore the natural 
energy structure of the problem. In [1], [9] and [10], Y. Guo and I. Tice introduced the geometric 
energy framework and prove both the local and global wellposedness in small data. The main 
idea in their proof is that instead of considering perturbation of the parabolic equation, we 
restart to prove the regularity under time-dependent basis. In our paper, we will follow this 
path and employ a similar argument. 

For the inviscid case, traditionally, we replace the no-slip condition with no penetration 
condition u-v = on It is often assumed that the initial fluid is irrotational and this curl-free 
condition will be conserved in the later time. This allows to reformulate the problem to one only 
on the free surface. The local wellposedness in this framework was proved by Wu [12, 13] and 
Lannes [14]. Local wellposedness without irrotational assumption was proved by Zhang-Zhang 
[15], Christodoulou-Lindblad [16], Lindblad [17], Coutand-Shkollcr [18] and Shatah-Zheng [19]. 
In the viscous case, vorticity will be naturally introduced at the free surface, so the surface 
formulation will not work. 

1.3. Geometric Formulation. In order to work in a fixed domain, we want to flatten the free 
surface via a coordinate transformation. Beale introduced a flatten transform in [2] and we will 
use a slightly modified version. We define a fixed domain 



for which we will write the coordinate x G fi. In this slab, we take E : {x 3 = 0} as the 
upper boundary and Eft : {X3 = —bo} as the lower boundary. Simply denote x' = {x\,x 2 ) and 
x = (xi, X2, xs). Then we define 



where V e 7] is defined as (A.3.1) when E = R 2 and as (A.4.1) when E = (L]T) x (L 2 T). 
Consider the transform from f2 to f2(i): 



(1.3.1) 



f2 = {x G E x R | - b < x 3 < 0} 



(1.3.2) 



rf = V e rj = The parameterized harmonic extension of r/ 



(1.3.3) 



$ e : {x 1 ,x 2 ,x 3 ) ^ {x 1 ,x 2 ,—x 3 + rf{\ + — )) = (yi, y 2 , y 3 ) 



This transform maps Q into Q(t) and its Jacobian matrix 



(1.3.4) 




1 






1 








and the transform matrix 



(1.3.5) 
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where 



A e = x^b/bo + d x ffb 
B e = x 3 d 2 b/b + d 2 fj e b 
(1-3.6) ,V = b/b + fj e /b + d 3 ffb 

K f = 1/J e 
b=l + x 3 /b 

By our assumption on initial data, there always exists 5 > 0, such that t]q(x') + b(x') > 2bo5 > 
for arbitrary x' G S. Based on lemma (A. 6) and lemma (A. 9), if we take e < 5 2 /(AC 2 ||?7o||#5/2), 
we have 



(1.3.7) 



< 5/2 



where tjq denotes if at t = 0. Furthermore, by fundamental theorem of calculus, we can estimate 
for any x 3 £ [—bo, 0] and 

(1.3.8) \fj e {x',x 3 )-fj e (x',0)\ < / \d 3 fj e (x',z)\d Z < b \\d 3 fj e \\ L oo <b S/2 

J—bn 



bo 

Let Jq denote J e at t = 0. Since f/o( x '>0) = Vo( x> ), w e have 

«/o = b/bo + Vo/b + d 3 fj e b=(b + rio)/b + (fj e 
(1.3.9) > 25-(5/2-V 2 = ^>0 



%)/^o + d 3 r] e b 



This means for given initial data i]o, we can always find < e < 5 2 /(4C 2 ||?7o||^5/2) such that 
Jq > 5 > 0, then this nonzero Jacobi implies the transform <& e makes sense at t = 0. If we 
further assume r/o £ H 7 ^ 2 (T,), we can conclude ^ is a C 1 diffeomorphism from f2 to Oo- 
In the following, we will always choose e = <5 2 /(4C 2 H7/0II//5/2) and for simplicity, we just write 
fj instead of rf , while the same fashion applies to A, <!>, A, B, J and K. For our special choice 
of e satisfying above condition, it is guaranteed that Jo > S > and $ is a C 1 diffeomorphism 
at t = 0. 

Define some transformed operators as follows. 

(V A f)i = Aijdjf 
V A - g = Aijdjgi 
A A f = • V A f 
N = (—dirj, -d 2 r], 1) 
{p A u)ij = Aikd k Uj + Ajkd k Ui 
S A {p,u) = pi — 3 A u 

where the summation should be understood in the Einstein convention. If we extend the diver- 
gence V_4- to act on symmetric tensor in the natural way, then a straightforward computation 
reveals that • S A (p, u) = V A p — A. A u for vector fields satisfying • u = 0. 
In our new coordinate, the original equation system (1.1.3) becomes 



(1.3.10) 



(1.3.11) 



dtu — dtr]bKd 3 u + u ■ V ' A u 
V4 • u = 
S A (p,u)M = f)N 
u = 

u(x, 0) = uq(x) 



A A u + V AP = 



m 
in 
on 
on 
in 



n 
n 

s 

n 



d t i] + uidi-q + u 2 <9 2 r/ = ^3 on S 

r)(x', 0) = i]o(x') on S 

where we can split the system into a Naiver-Stokes equation and a transport equation. 
Since A is determined by r] through the transform, so all the quantities above is related to f], 
i.e. the geometric structure of the free surface. This is the central idea of our proof. It is 
noticeable that in proving local wellposedness of above equation system, we must verify <l?(t) is 
a C 1 diffeomorphism for any t G [0, T], where the theorem holds. 
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1.4. Main Theorem. In this paper, we will prove the local wellposedness for higher order 
regularity. 

Theorem 1.1. Let N > 3 be an integer. Assume the initial data rjo + b > 2&o<5 > 0. Suppose 

1 1 2 //j 

i/iai wo anc ^ Vo satisfy the estimate \\uo\\ H 2n + ||?7o||#2iv+i/2( S ) < oo as well as the N compatible 
condition. Then there exists < To < 1 such that for < T < Tq, there exists a unique solution 
(u,p,r)) to system (1.3.11) on the interval [0, T] that achieves the initial data. Furthermore, the 
solution obeys the estimate 




< C(Q ,S)P(\\u \\ H 2N + |MIh2JV+1/2( S) ) 

where C(Qq, 5) is a positive constant depending on the initial domain fig and separation quantity 
5 and P(-) is a single variable polynomial. The solution is unique among functions that achieves 
the initial data and for which the left hand side of the estimate is finite. Moreover, rj is such 
that the mapping <£(i) defined by (1.3.3) is a C 2N ~ 2 diffeomorphism for each t G [0, T\. 

Remark 1.2. In above theorem, H k denotes the usual Sobolev space in Q and denotes 
the usual Sobolev space on S, while we will state the definition of \\-\\ x * later. 

Remark 1.3. Because the map $>(t) is a C 2N ~ 2 diffeomorphism, then we may change variable 
to y £ Q(t) to produce solution of (1.1.3). 

Remark 1.4. Since our paper is a natural extension of results in [1] from small data to arbitrary 
data and we employ the same framework, then many lemmas and theorems are quite similar 
except for some minor corrections. Hence, we will not always give the detailed proof of these 
results. Instead, we will refer to the corresponding results in [1] and merely point out the 
differences when necessary. 

1.5. Convention and Terminology. We now mention some of the definition, bits of notation 
and conventions we will use throughout the paper. 

(1) We will employ Einstein summation convention to sum up repeated indices for vector 
and tensor operations. 

(2) Throughout the paper C > will denote a constant only depend on the parameter of 
the problem, N and Q, but does not depend on the data. They are referred as universal 
and can change from one inequality to another one. When we write C(z), it means a 
certain positive constant depending on quantity z. 

There are two exceptions to above rules. The first one is that in the elliptic estimates and 
Korn's inequality, there are constants depending on the initial domain Qq. Although this 
should be understood as depending on the initial free surface 770, since its dependent 
relation is given implicitly and cannot be simplified further, we will also call them 
universal. The second one is that we will also call the constant C(S) universal, though 
5 is determined by initial domain Qq. To note that apart from these, all the other 
constants related to initial data Qq, uq and t]q should be specified in detail. 

(3) We will employ notation a < b to denote a < Cb, where C is a universal constant as 
defined above. 
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(4) We will write P(z±, . . . , z n ) to denote the polynomial with respect to finite arguments 
Zi for i = 1, . . . , n. The detailed properties of the polynomial will be specified in the 
context. This notation is used to denote some very complicated polynomial expressions, 
however whose details we do not really care. This kind of polynomial may change from 
line to line. 

(5) We always write Df to denote the horizontal derivative of / and V/ for the full deriv- 
ative. 

(6) For convenience, we will typically write H° = I? , except for notation L 2 ([0,T]; H k ). 
We write H k (Q,) with k > and H S (T,) with s 6 i! for standard Sobolev space. Same 
style of notation also holds for L 2 ([0,T];H k ) and L°°([0, T};H k ). When we write ||-|| Hfc , 
this always means the Sobolev norm in Q, otherwise, we will point out the exact space 
it stands for, e.g. IHIijfc(s)- ^ similar fashion is adopted for ||-|| L2 ^fc and ||-|lL°°_f/ fc - 

1.6. Structure of This Paper. Our proof of the main theorem is based on an iteration 
argument for nonlinear terms. Hence, in Section 2, we will first prove the wellposedness to 
higher regularity for the linear equation 

dtu — A^n + V ' A u = F in fl 

V.4 ■ u = in Q 

(pi - B A u)M = H on S 

u = on Sb 

u(t = 0) = no in Q 

In Section 3, we prove the wellposedness for transport equation 

dtV + u\d\r\ + n 2 c>2f? = n 3 on S 
i](t = 0) = i]o on S 

and also the estimate for forcing terms in nonlinear problems. Finally, in Section 4, we will 
employ an iteration argument to complete the proof. 

Throughout the paper, we always assume N > 3 is an integer. We consider both infinite case 
and periodic case simultaneously and won't differentiate between them unless necessary. 



(1.6.1) 



(1.6.2) 



2. Linear Navier-Stokes Equation 

2.1. Introduction. In this section, we will concentrate on proving the wellposedness of the 
linear Navier Stokes equation 



d t u 


- A A u + V^n = F 


in 


n 


V.4 


■ u = 


in 


n 


(pi 


- B A u)M = H 


on 


s 


u = 





on 


s b 


u(t 


= 0) = n 


in 


n 



Throughout this section, we will assume for any t in the finite time interval [0, T] which we will 
specify later, the condition J(t) > 5/2 > holds, where J(t) is the Jacobi of transform <&(£). 

2.2. Preliminaries. 

2.2.1. Transform Estimates. In order to study the linear problem in the slab domain, we will 
employ the idea that transforming the variable-coefficient problem into a constant-coefficient 
problem through diffeomorphism. So before estimating, we need to confirm the mapping $ is 
an isomorphism from to £l(t) and determine the relation of corresponding norms between 
these two spaces. We borrow the idea from lemma 2.3.1 of [1]. 

Lemma 2.1. Let *S> : Q — > Q' be a C 1 diffeomorphism satisfying ^ G H^ 1 , the Jacobi J = 
det(V^) > 5 > a.e. in n and W - / G H k (Q) for an integer k > 3. If v € H m (Q r ), then 
v o * g H m (n) form = 0,1,..., k + 1, and 

(2-2.1) ||n o ^\\ H m (n) < Ci|H| Hm{f y) 



6 
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Gi = (l + ^)(l + ||V*-7|| Hfc(n) ) + 



Similarly, for u G H m (n), wo^e /or m = 0, 1, . . . , jfe + 1 and 

(2-2.3) ||u o ^~ 1 || H m (n , ) < C2||«||H"*(n) 



(2.2.4) 



C 2 



l + ^)( l + ||Vvr/-/| 



\ 2fc+2 



Let £' = \I/(E) denote the upper boundary of SV . If v G #™-V2(vy) /or m = 1, . . . , A; — 1, t/ien 

(2.2.5) O || H m-l/2( E ) < C 3 1^11^-1/2(2,) 

where 

(2.2.6) C 3 = C 1 C(||V*-/|| Hfc(n) ) 

//we ,H" m - 1 /2( S ) form = l,...,k-l, thenuo^~ l G fl™" 1 ^!/) anrf 

(2.2.7) ||u o ^~ 1 ||^ m _ 1/2( - E ,- ) < C , 4||u|| ffm -i/2( E ) 
where 

(2.2.8) C 4 = C 2 C(||V*-/|| Hfc(n) ) 

/n all above, C(\\V^ — I\\jjk^) is a constant depending on \\V^ — I\\jjk^y 

Proof. Let x denote the coordinate in f2 and y in CI'. A direct computation shows that, for the 
case m = 0, if we have v G L 2 (Q'), then 



\v o ty\ 



L 2 (Q) 



[ \voV(x)\ 2 dx= [ Ivfjixy'dyKpix)- 1 ]] I b| 2 dy 

For m = 1, G L 2 (S1'), then 
(2.2.9) 



— (t, o = — o *( x )_i( s ) 



<9x 



Hence, we have the estimate 



d_ 

dx. 



(v o *) 



L 2 (Q) 



<9f 



9*, 



dxi 



< 1 + 



< 1 + 



<9v 



1 + 



<9f 



L 2 (f2') 



9*, 



(s) 



L°°(f2) 



L 2 (f2') 



< 



r 


9w 







L°°(n) 



L 2 (a') 



< (l + J)(l + l|V*-/|| H *(n)) ll«llHi ( no 



Similarly, it is easy to show the result when m = 2, i.e. if d\v G H (O'), then 

2 / i \ / \ 4 



5 , TX 



if!(n) 



< (l + J) (l + II V* - /|| H * ( n) ) \M 2 ma ,) 
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For m = 3, . . . , k + 1, inductively, we can assume the statement is valid for m < tuq and try to 
justify it in ttiq + 1 case. 



d dv d^ ' dv dv 

(2.2.10) — (v o = — o ^( x )_2( x ) = — o ^r(x) + — o ^(x) 



dxi 



(x) - h 



Induction hypothesis tells us if v G H m ° +1 , then we have 

<% 

Since we have the multiplicative embedding ff m o x H k if mo for itlq > 2 and k > 3, we 
deduce that 

i«o$) G ff m ° 
ox, 

which implies oo$£ /-fm +i_ Moreover, we have the estimate 



dx. 



(v o *) 



_ff m o(a) 



< 



< 



< 



dv 



+ 



ff m o (SY) 



dv 



1 + IIV* - f| 



1 



^ T / \ 



<9x; 

2 



(x) - f, 



H m o (CI) 



1 + T l + ||W-f| 



2m +2 



\ v \\H rn o+ 1 (n>) 



\ m-o+l 



lij^o+^n') 



which implies that 

Ik *ll/r»o+i(n) < (* + ]) ( x + II W " 7 Hff*(n) ) 

Therefore, (2.2.1) holds. A similar argument justifies (2.2.3), utilizing the fact V\I/ _1 (y) = 
(V*)- 1 o*-!(y). 

Since G H^ 1 , we have £' is locally the graph of a C k ~ 1 ' 1 / 2 function. Based on [6], there 
exists an extension operator E : if m ~ 1 / 2 (Tl') — > H m (£l') for m = 1, . . . , k — 1 with the norm of 
the operator depending on C(||W - 7|| H * (n) ). For v G ff m - 1 / 2 (s / ), let V = f5t> G H m {n'). 
Furthermore, we have V o ^ G H m (Q) and i;ovI/ = V r ovt'| s g jjm-i/2 Moreover, 

(2-2.11) ||vO^|| Hm _ 1/2(s) < ||F0^|| Hm(n) < Cl||F|| H m( Q ,) < C 3 ||w|| Hm -i/2 (E ,) 

which is (2.2.5). Similarly, (2.2.7) follows. □ 

Remark 2.2. Based on our assumptions on fj and b, it is easy to see $ defined in 1.3.3 is a 
C 1 diffeomorphism satisfying the hypothesis in lemma 2.1. Moreover, the universal bounding 
constant C now can be taken in a succinct form 



(2.2.12) 



c = c(lM| Hfc +i/2 (s) ) 



to replace Cj in the estimate. If we only consider the estimate in the domain f2 and need to 
specify the constant, we can take the explicit form 



(2.2.13) 

to replace C\ and Ci- 



c = {i + H 



tffc+l/2( S ), 



\2fc+2 



2.2.2. Functional Spaces. Now we introduce several functional spaces. We write H k (Q) and 
if fe (S) for usual Sobolev space of either scalar or vector functions. Define 



(2.2.14) 



W(t) = {u(t)eH 1 (n): u(t)| Sfc =0} 

X(t) = {u(t) G ff 1 ^) : «(t)| E6 = 0, V4 (t) • u(t) = 0} 



LEI WU 



Moreover, let L 2 ([0, T};H k (fl)) and L 2 ([0, T];H h (T,)) denotes the usual time-involved Sobolev 
space. Define 



(2.2.15) 



W = {u€L 2 ([0,T};H 1 (Q)) : u| Sfc = 0} 

X = {ueL 2 ([0,T};H\n)) : u| Sfc = 0, • u = 0} 



It is easy to see W, X, W and X are all Hilbert spaces with the inner product defined exactly 
the same as in H l {Q) or L 2 ([0,T]; H 1 ^)). Hence, in the subsequence, we will take the norms 
of these spaces the same as || • \\ H i or || • ||^2^i. 

Furthermore, we define the functional space for zero upper boundary. 

(2.2.16) V{t) = {u{t) e H 1 ^) : u(t)| E = 0} 

Also we will need an orthogonal decomposition H°(Q) = Y(t) © Y- 1 ^), where 

(2.2.17) Y ± (t) = {V A(t) <p(t) : <p{t) G V(t)} 

In our use of these spaces, we will often drop the (t) when there is no potential of confusion. 
Finally we will define the regular divergence- free space in H l (Q). 

( X = {ueH 1 (Q) : n| Eb = 0,V-n = 0} 

1 j X = {u€L 2 ([0,T];H 1 (n)) : n| Efc = 0, V • u = 0} 

Sometimes, we will consider the functional space with the similar properties as defined above 
in transformed space il'. If that is the case, we will employ the notation W(£l') or X(£l') to 
specify the space they live in. The same fashion can be applied to other functional spaces. 
We will use (■, -)-^o to denote the normalized inner product in H°, i.e. 

(2.2.19) (u,v) n o = / Ju-v 

Jn 

and (■, -}c 2 h to denote the normalized inner product in L 2 ([0, T]; H°), i.e. 

(2.2.20) (u,v) C 2 H o= [ [ Ju-v 

Jo Jn 

Also we will use (•, -)w and (■, -)w to denote the regular inner product in H l if both arguments 
are in If or W. Moreover, we define the regular inner product on £ as follows. 



(2.2.21) («,w)wo (E) 



u ■ v 



E 



(2.2.22) (u,v)& H op) =11 

JO JY. 



u ■ V 

s 



The following result gives a version of Korn's type inequality for initial domain fio- Here we 
employ Beale's idea in [2]. 

Lemma 2.3. Suppose Qq is the initial domain and ryo G H 5 ^ 2 (T,). Then 
(2-2-23) \H 2 H H n ) Z W^\\ 2 H o m 

for any v G // 1 ($7o) and v = on {ys = —b}. 

Proof. In the periodic case, S is bounded and v\-£ b = 0, so naturally Korn's inequality is valid 
(see proof of lemma 2.7 in [2]). Then we consider the infinite case. First we prove the decaying 
property of initial surface 770, i.e. for \a\ < 1, |9 Q r/o| — > as \x\ — > 00 in the slab. For horizontal 
derivative d a , 



d a f} = [ {2Ttii) a e ex ^e 2 ™'-%{Z)&t 
Jr? 
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where rjo is the Fourier transform of rjo in R 2 . Then 
/ (2m0 a e^U0dt< [ \£\ a e^\UO\a + \t\ 5/A )( 

Vi + |£| 5/ / 

The last inequality is valid since X3 < 0. Hence, (2m£) a e eX3 ^\rio(£) G L l (R 2 ). A similar proof 
can justify the result if d a = 83. Since \x\ — > 00 naturally implies \x'\ — > 00 in the slab, the 
Riemann-Lebesgue lemma implies |9 a r/o| — > as \x\ — > 00. 

We construct a mapping a = $(0) as defined in (1.3.3), which maps f2 = {R 3 : — feo < £3 < 0} 
to Qo. Denote a(x) = x + a(x). The above decaying property and our assumption on b leads 
to d a a(x) — > as \x\ — > 00 for |a| < 1. 
We partition the slab f2 into cubes 

Qj,k = {x : j < xi < j + 1, k < x 2 < k + 1, -6 < £3 < 0} 

In each cube, we have Korn's inequality 

IMItfl(Q) < C , (||Bt>||^ 0( Q ) + ||v||^o (Q) ) 

Employing the compactness argument as Beale did, under the condition v = on {X3 = — 60}, 
we can strengthen this result to 

IMIjyi(Q) < C \\Bv\\ 2 H0{Q) 

This argument relies on the fact that for such v 

\\v\\ 2 H o {Q) < C \\Bv\\ 2 H o {Q) 

Now suppose D C £l is the image of union of such cubes contained in {\x\ > R} for some R. 
Applying last estimate to v = uo a and then transform to Slo, we have 

IMItfi(D) - C 1 1 ^ 1 1^0(0) + 6 IMIfri(D) 

where e < 1 when R is large enough based on the decaying property. Then we use Korn's 
inequality in the bounded domain Q$ — D and combine with above estimate to obtain the 
estimate required. □ 

Remark 2.4. In the final step of the proof, the bounding constant in Korn's inequality depends 
on the domain Qq. However, based on our notation rule, we still use < to denote it. 

Next, we need to show the equivalence of certain quantities in Q. 

Lemma 2.5. There exists a < eo < 1, such that if 1 1 77 — ?7o || < eo, then the following 
relation 

(2.2.24) |H|^„ < J J \u\ 2 < (1 + |MI;j5/2 (s) ) IMI^o 



3 II ||2 



(2-2.25) I ^ ||n||^ < / J |D^| 2 < (1 + ||%|| H 5/ a(E) ; 

holds for all u £ W. 

Proof. (2.2.24) is just a natural corollary of relation 5 < ||J|| L oo < (1 + ||V^|| ioo ) < (1 + 
ll r ?ll_H'5/2) < (1 + 1 1 ??o 1 1 ) based on our assumption and Sobolev embedding. 
To derive (2.2.25), notice that 

(2.2.26) f J\O a u\ 2 = f J\B Ao u\ 2 + [ J(0 A u + A() u) : (B A u-D Ao u) 
Ju Jci Jo. 
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For the first term on the right-handed side, based on the integral substitution and Korn's 
inequality we proved above, we have 

/ J\B Ao u\ 2 dx > ) / J |IBU o u| 2 dx= \ / |Bi;| 2 dy 

Jn 1 + 11% 1 1 #5/2^ Jn i + 1 1 % 1 1 #5/2^ Jn 



1 



> IhillL ,„ \ > 



l + ho||# 5 / 2(s) l|V ^ 1( ^-(l + || % ||# 5 / 2(s) ) 3 lMHl 



For the second term on the right-handed side, using Korn's inequality in slab fi, we can naturally 
estimate 

[ J(D A u + B Ao u) : (B A u-0 A() u) < \\J\\ LOO \\A + A \\ L °o \\A- A \\ L oo [ |Vu| 2 
Jn Jn 

< e (l + |M|#5/ 2(S) ) 3 jf |Bn| 2 < e (l + ||%||#5/ 



2(E)) 3 ll n 



For given initial data, we can always take eo sufficiently small to absorb the second term into 
the first one. Then we have 

/ J|B^| 2 > r I, ,} 73 ||u||^i 

Jn \± + 1 1 % 1 1 #5/2(2) J 

The first part of (2.2.25) is verified. For the second part, notice that 

/ 2 / 2 2 2 2 

/ J\^au\ < (1 + IMI#5/ 2 ) / \B A u\ < (1 + ||77o|| H 5/2)max{lJ|AA'|| LOO , ||Biir|| LOO , ||if|| LOO } ||u|| 
Jn Jn 

Since it is easy to see 

max{l, \\AK\\ 2 LOO , \\BK\\ 2 Loa , ll^lllco} < 1 + (1 + llV^Hico) H^Hloo < (1 + ||%||#5/2 (s) ) 2 
then the second part of (2.2.25) follows. 

□ 

Remark 2.6. Throughout this section, we can always assume the restriction of rj depending on 
eo is justified, and finally we will verify this condition for t G [0, T] in the nonlinear part. 

Now we present a lemma about the differentiability of norms in time-dependent space. 

Lemma 2.7. Suppose that u G W and dtu G W*. XTien i/ie mapping t — > \\u(t)\\^jo is absolute 
continuous, and 

(2.2.27) ^ ||n(t)||^o = 2(^u(t), u(t)> + \u(t)\ 2 8 t J(t) 
Moreover, u G C°([0,T]; H°(Q)). If v € W and d t v G W* as well, we have 

(2.2.28) A /" „(t) . „( t ) = (5 t u(t), v(t)) + {d t v(t), u(t)) + / u(t) • v(t)d t J(t) 

Proof This is exactly the same result as lemma 2.4 in [1], so we omit the proof here. □ 

Next we show the estimate for H~ 1 / 2 boundary functions. 

Lemma 2.8. // v G H° and X7 A • v G H°, then v • M G H' 1 / 2 ^), v ■ u £ H' 1 / 2 ^) (1/ is the 
outer normal vector on ,) and satisfies the estimate 



(2.2.29) 



|« • A/"|| ff -l/2 (S) + \\V ■ v\\ H -l/2p b ) ^ (! + ll"ll#5/2) 2 ^ HvllflO + ||V^ • V\\ H ^ 
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Proof. We only need to prove the result on S; the result on £5 can be derived in s similar 
fashion considering JAe^ = v on £5. 

Let <p G .ff 1//2 (£) be a scalar function and let <fi G V be a bounded extension. We have 
/ (j)vN = / JAijVi^ej ■ e 3 ) = / V4 • (v<j>)J = / <^V^ -vJ + v V A 4>J 

v A 4> 



< (i + hWII^ /2 )( <p iy0 ||v >A -«|| fl t.+ 



\H° 



H° 



< 



(l + \\v(t)\\ H ^) (\Mh° + \\Va-v\\ ho) 11^11^1/2(2) 



Since <fi is arbitrary, our result naturally follows. □ 

Remark 2.9. Recall the space y(t) C H°. It can be shown that if v G y(t), then V A • v = 
in the weak sense, such that lemma 2.8 implies that v ■ M G ff _1 / 2 (£) and v ■ v G i? -1 / 2 (£;,). 
Moreover, since the elements of y(t) are orthogonal to each V A (p f or 4> £ ^"(*)> roe Z" 1 ^ 
■u • v = on £&. 

We want to connect the divergence-free space and divergence-A-free space, so we define 

/ K 

(2.2.30) M = M(t) = KV<S> = K 

\ ART Bif 1 

Note that M is invertible and M _1 = JA T . Since dj(JAij) = for i = 1,2,3, we have the 
following relation 

p = V A - v o Jp = JV A ■ v = JAijdjVi = djiJAijVi) = dj(JA T v)j = d j (M~ 1 v) j = V • (M^v) 

Then if V_4 • v = 0, we have V • (M~ l v ) = 0. M induces a linear operator Mt : u — )■ .Mt('u) = 
M(t)u. It has the following property. 

Lemma 2.10. For each t G [0, T], A; = 0, 1, 2, .Mi is a bounded linear isomorphism from H k (Q) 
to H k (£l) and from Xq to X. The bounding constant is given by (1 + ||?7(i)||#-9/2( S )) 2 . If we 
further define M. defined by Mu(t) = M. t u(t), then it is a bounded linear isomorphism from 
L 2 ([Q,T]]H k {9)) to L 2 ([0,T];H k ({})) and from Xq to X. The bounding constant is given by 
(1 + sup < t < T ||r/(t)|| H 9/2 (s) ) 2 . 

Proof. It is easy to see for each t G [0, T], 

\\M t u\\ Hk < \\M t \\ C 2 \\u\\ Hk < (1 + ||r?(t)||^9/2 (s) ) 2 \\u\\ Hk 

for k = 0, 1, 2, which implies Ai t is a bounded operator from H k to Since .Mt is invertible, 
we can estimate < (1 + ||??(i)||//9/2(2)) 2 IMIiffc- Hence Mt is an isomorphism be- 

tween H k . Also above analysis implies Mt maps divergence- free function to divergence- A- free 
function. Then it is also an isomorphism from Xq to X(t). 

A similar argument can justify the case for L 2 ([0, T];H k ) and X. □ 

2.2.3. Pressure as a Lagrangian Multiplier. It is well-known that in usual Navier-Stokes equa- 
tion, pressure can be taken as a Lagrangian multiplier. In our new settings now, proposition 
2.2.9 in [1] gives construction of pressure from transformed equation 1.3.11, which is valid for 
small free surface. So our result in arbitrary initial data is not completely ready, which we will 
present here. 

For p G H°, we define the functional St G W* by St(v) = {p,V A • v)yo. By the Riesz represen- 
tation theorem, there exists a unique Qt{p) G W such that St(v) = {Qt(p),v)w for all v G W. 
This defines a linear operator Q t : H° — > W, which is bounded since we may take v = Qt(p) to 
see 

\\Qt(p)\\w = S t (v) = ( P ,V A ■ v) H o < (1 + 11^)11^5/2(2)) \\ P \\ H0 \\V A ■ v\\ H0 



< 



(1 + ll^)ll^/ 2(s) ) 3 \\p\\ H o \\v\\ w = (1 + ||»7(t)|| H 5/2 (E) ) 3 HpIIho \\Qt(p)\\ w 
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so we have ||Qt(p)||w < (1 + \\v(t)\\H 5 / 2 (z)) 3 \\p\\h°- Similarly, for S G W*, we may also define a 
bounded linear operator Q : L 2 ([0,T]; H°) — > W via the relation {p,^A' v )c 2 v. = (Q{p)i v )w = 
S(v) for all neW. Similar argument shows ||Q(p)|| w < (1 + sup < t < T ||r?(*)ll.ff5/2(E)) 3 l|p|lz, 2 .ff°- 

Lemma 2.11. Ze£ p G i/ten i/iere exists a v £ W such that • v = p and \\v\\ w < 
(1 + ||^(i)||^5/2( S )) 6 ||p||fl-o- 7/ instead p G L 2 ([0, T]; if ), i/jen i/iere exists a f G W suc/i i/tai 

tf 5 /2(£)) 



VA-v=p for a.e. t and \\v\\ w < (1 + sup 0<t<T ||r7(t)||^ 5 /2 (2) ) 6 ||p|| L 2 H o 



Proof. In the proof of lemma 3.3 of [2], it is established that for any g G L 2 (f2), the problem 
V • u = q admits a solution u G W such that < IMI//0- A simple modification of this 

proof in infinite case can be applied to periodic case. Let define q = Jp, then 

hWlo = J q \q\ 2 = J a \p\ 2 J 2 < || Jlli. ||p|||o < (1 + ||r?(t)|| ff6 / 2(s) ) 2 |b||^o 
Hence, we know v = M(t)u G W satisfies V^w = p and 

\Mw ~ ( 1 + ll»7(*)llif9/2(s)) 4 |l u llw~( 1 + ll ? 7(*)lli?9/ 2 (s)) 4 |kll^o 

< (i + II^)II^/ 2( e)) 6 Nho 

A similar argument may justify the case for W. □ 

With this lemma in hand, we can show the range of operator Q t and Q is closed in W and 
W. 

Lemma 2.12. R(Qt) is closed in W , and R(Q) is closed in W. 

Proof. For p G if , let v G be the solution of V^ ■ v = p provided by lemma 2.11. Then 

\\p\\h<> ~ (P^A • v)yp = {Qt(p),v)w < IIQt(p)llw \\ v \\w IIQt(p)ll W (l+ll»/(t)llfl»/2(E)) 6 Wflt. 
such that 

(1 + [1^)11^^)3 ll&GOHw £ llPllflf ^ C 1 + lk(t)llflt>/2 (E )) 6 \\Qt(p)\\ w 

Hence R(Qt) is closed in W. A similar analysis shows that i?(Q) is closed in W. □ 

Now we can perform a decomposition of W and W. 

Lemma 2.13. We have that 1^ = 1® R{Qt), i-e. X 1 - = R(Q t ). Also, W = X R(Q), i.e. 
X 1 - = R(Q). 

Proof. By lemma 2.12, R(Qt) is closed subspace of W, so it suffices to show R(Q t ) ± = A. Let 
v G R(Qt) ± , then for all p G if , we have 

(p,V A ■ v) H o = (Qt{p),v) w = 

and hence ■ v = 0, which implies R{Qt) 1 ' C A. On the other hand, suppose v G A, then 
V^4 ■ w = implies 

(Qt(p),^)w = (p, • v) H o = o 

for all p G if . Hence V G f^Qt) 1 " and we see X ^ R(Q t ) ± . So we finish the proof. A similar 
argument can show W = X ® R{Q). □ 

Proposition 2.14. If X £ W* such that X(v) = for all v G A, then there exists a unique 
p(t) G if such that 

(2.2.31) (p(t), V4 • v) w o = A(v) /or all v £ W 
and \\p(t)\\ H o < (1 + ||r/(i)||ff9/ 2(s) ) 6 ||A||w-. 

7/ A G W* suc/i i/iai A(u) = /or aZZ v £ X , then there exists a unique p £ f 2 ([0, T]; if ) such 
that 

(2.2.32) (p, V.4 • v) C 2 H o = A(v) for all v £ W 
an d \\p\\ L 2 H o < (1 + sup < t < T ||r/(0|| H g/ 2(s) ) 6 ||A||w» 
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Proof. If X(v) = for all v £ X, then the Riesz representation theorem yields the existence of 
a unique u G X 1 - such that X(v) = (u,v)w for ah v G W. By lemma 2.13, u = Q(p) for some 
p G H°. Then A(v) = {Qt(p),v)w = {p(t),V A ■ v) n o for all v G W. 

As for the estimate, by lemma 2.11, we may find v e W such that • v = p and \\v\\ w < 
(l + ll^)llfl«/2(E)) 6 WflO- Hence > 

IIpIIho < (p,V^-w) W o = A(v) < \\\\\w*(l + \\v(t)\\H*/Hs/ Mho 
So the desired estimate holds and a similar argument can show the result for A. □ 

2.3. Elliptic Estimates. In this section, we will study two types of elliptic problems, which will 
be employed in deriving the linear estimates. For both equations, we will present wellposedness 
theorems to the higher regularity. 



2.3.1. A-Stokes Equation. Let us consider the stationary Navier-Stokes problem. 
(2.3.1) 



V.4 • S A (p, u) = F in Q 

V A • u = G in $7 

Sa(p, u)M = H on T 

u = on Tii 



Since this problem is stationary, we will temporarily ignore the time dependence of n,A, etc. 
Before discussing the higher regularity result for this equation, we should first define the weak 
formulation. Our method is quite standard; we will introduce p after first solving a pressureless 
problem. Suppose F G W*, G G H° and H G H' 1 / 2 ^). We say (u,p) G W x H° is a weak 
solution to 2.3.1 if • u = G a.e. in ft, and 

(2.3.2) ^AU, n A v) W) - (p, V A ■ v) n o = (F, v) n o - (H, v) H0(s) 
for all v G W. 

Lemma 2.15. Suppose F G W* , G G H° and H G i^~ 1 / 2 (S), then there exists a unique weak 
solution (u,p) £W x H° to 2.3.1. 

Proof. By lemma 2.11, there exists a u G W such that • u = F 2 . Naturally we can switch 
the unknowns to w = u — u such that in the weak formulation w is such that ■ w = and 
satisfies 

(2.3.3) ^(B A w,B A v) H o - (p,V A ■ v) n o = ~<Pau,B a v) u o + (F, v) H o - (H,v) H o { ^ 
for all v G W. 

To solve this problem, we may restrict our test function to v G X such that the pressure term 
vanishes. A direct application of Riesz representation theorem to the Hilbert space whose inner 
product is defined as (u,v) = (D.4U, B^f ) w o provides a unique w G X such that 

(2.3.4) ^<p A w,B A v) H o = -^(B A u,D A v) n o + (F,v) H o - {H,v) H op) 
for all v G X. 

In order to introduce the pressure, we can define A G W* as the difference of the left and right 
hand sides in (2.3.4). So X(v) = for all v£l. Then by proposition 2.14, there exists a unique 
p G H° satisfying (p, V A • v)yo = X(v) for all v G W, which is equivalent to (2.3.3). □ 

The regularity gain available for solution to (2.3.1) is limited by the regularity of the coeffi- 
cients of the operator A_4, V.4 and V^-, and hence by the regularity of n. In the next lemma, 
we will present some preliminary elliptic estimates. 

Lemma 2.16. Suppose that n G H k+1 / 2 (T) for k > 3 such that the map $ defined in (1.3.3) is 
a C l diffeomorphism of SI to Q' = If F G H°{n), G G H 1 ^) and H G H 1 / 2 ^), then the 
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equation (2.3.1) admits a unique strong solution (u,p) G H 2 x H 1 , i.e. (u,p) satisfies (2.3.1) 
in strong sense. Moreover, for r = 2, . . . , k — 1 we have the estimate 



(2.3.5) 



\ U \\m + lbllff-1 < C(??)f ll^ ? ll//r-2 + + \\H\\ H r-3/2^ 



whenever the right hand side is finite, where C(n) is a constant depending on \\n\\ H k+i/2^y 

Proof. This lemma is exactly the same as lemma 3.6 in [1], so we omit the proof here. □ 

Notice that the estimate (2.3.5) can only go up to k — 1 order, which does not fully satisfy 
our requirement. Hence, in the following we will employ approximating argument to improve 
this estimate. For clarity, we divide it into two steps. In the next lemma, we first prove that 
the constant can actually only depend on the initial free surface. 

Lemma 2.17. Let k > 6 be an integer and suppose that n G tf fe+1 /2(£) and % G # fe+1/2 (£). 
Then there exists eo > such that if \\n — rjo\\ H k-3/2r^ < eo, solution to 2.3.1 satisfies 

(2.3.6) \\u\\ H r + IbHifr-l < C(r]o)( \\F\\ H r-2 + \\G\\ H r-l + \\H\\ H r-3/2 



2 (S) 

for r = 2, . . . , k — 1, whenever the right hand side is finite, where C(t]q) is a constant depending 

071 11% 1 1^ + 1/2(2)- 

Proof. Based on lemma 2.16, we have the estimate 

(2.3.7) \\u\\ H r + ||p||#r-l ^ C(q) ( \\F\\ Hr -2 + \\G\\ H r- 1 + \\H\\ Hr -- i /2 



2 (2) 

for r = 2, . . . , k — 1, whenever the right hand side is finite, where C(rj) is a constant depending 
on |M| H fc+i/2 (E ). Define £ = r? - r? , then £ G H k+1 / 2 (E). 

Let us denote Aq and A/"o of quantities in terms of rjo. We rewrite the equation 2.3.1 as a 
perturbation of initial status 



(2.3.8) 



where 



{ V Ao -S Ao {p,u) = F + F° in n 

V Ao -u = G + G° in ft 

S Ao (p,u)M = H + H° on E 

u = on T,b 



F° = V Ao - A ■ S A (p, u) + V^ • S Ao - A (p, u) 

(2.3.9) G° = V Ao - A • u 

H° = S Aq (p, u)(M -M)+ S M - A (p, u) 

Suppose that ||£ \\ H k-m{T.) < 1, which implies ||£||^fc-3/2( S ) < ||£||H fc - 3 / 2 (E) < 1 f° r an y ^ > 1- A 
straightforward calculation reveals that 

H-^llflr-a < C(l + ||?7o||iffe+i/2 (E )) 4 ||£||ff*-3/2 (E) (\\u\\ H r + |b||ifr-l) 

(2.3.10) H^l^-i < C(l + |MU+i/2 (s) ) 2 ||eilH*-3/2 (E) (\\u\\ Hr ) 

\\H\\ H r-3/2^ < C(l + ||r/o|| H fc + l/2( S )) 2 11^11^-3/2(2) + INI^-l) 

for r = 2, . . . , k — 1 and C includes the term related to <5. 

Since the initial surface function 770 satisfies all the requirement of lemma 2.16, we arrive at the 

estimate that for r = 2, . . . , k — 1 

(2.3.11) 

h\\ H r + \\p\\ H r-i < CM ( \\F + F°|| Hr - 2 + \\G + G°\\ Hr _ 1 + \\H + H Q \\ Hr _, /2(J]) ) 
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where C(%) is a constant depending on ||?7o||#k+i/2(£)- Combining all above, we will have 
(2.3.12) 

\\ u \\h- + Ibllfr-i ~ 

C(V0) ^ \\F\\ H r-2 + \\G\\„r-l + ||#||ffr-3/2 (S) ^ + C{r] Q ){\ + 1 1 T] \ | ^fc+1/2 (S) ) 4 ||£||ffk-8/a (E) 

So if 

(2.3.13, K || ff ^ (E) < vcMil + M^* } 

we can absorb the extra term in right hand side into left hand side and get a succinct form 

(2.3.14) \H\ H r + IblljF-l < C(7] )(^ \\F\\ Hr -2 + ||G|| H r-l + || H || H r-3/2( S) 

forr = 2,...,/c-l. □ 

Note that the above lemma only concerns about regularity up to k — 1 and we actually need 
two more order. Then the next result allows us to achieve this with a bootstrapping argument. 

Proposition 2.18. Let k > 6 be an integer, and suppose that n £ H k+1 ^ 2 (T,) as well as 
r/o G H k+l l 2 {Yi) satisfying \\n — r]Q\\ H k+i/2(^) < e o- Then solution to 2.3.1 satisfies 



(2.3.15) \\u\\ H r + \\p\\ H r-i < C(rjo) ( ||-F|l H r-2 + \\G\\ Hr -i + \\H\\ Hr - 3 /2^ 

for r = 2, . . . , k + 1, whenever the right hand side is finite, where C(r/o) is a constant depending 

on ll r /o||/ffc+i/2( S )- 

Proof. If r < k— 1, then this is just the conclusion of lemma 2.17, so our main aim is to gain two 
more regularity here for r = k and r = k + 1. In the following, we first define an approximate 
sequence for n. In the case that £ = R 2 , we let p 6 Cq°(R 2 ) be such that supp(p) C -8(0,2) 
and p = 1 for 5(0, 1). For m £ N, define n m by r) m (£) = p(£,/m)fj(£) where ~ denotes the Fourier 
transform. For each m, n m £ iiP(X) for arbitrary j > and also n m — > n in ff fc+1 / 2 (E) as 
m — > co. In the periodic case, we define ff 1 by throwing away the higher frequencies: f} m (n) = 
for |n| > m. Then n m has the same convergence property as above. Let A m and N m be defined 
in terms of n m . 

Consider the problem (2.3.1) with A and M replaced by A m and M m . Since rf 1 G H k+5 / 2 , we 
can apply lemma (2.16) to deduce the existence of (u m ,p m ) that solves 

' V Am ■ S A m (p m , u m ) = F in n 
(2316) V A -u™ = G in n 

u m = on £ 6 

and such that 

(2.3.17) ||«™|| Hr + \\p m \\ H ^ < C(||n m ||^ + 5/ 2(E) ) ( ||F|| Hr _ a + \\G\\ Hr -i + ||# ||^-3/ 2(E) ) 

for r = 2, . . . , k + 1. We can rewrite above equation in the following shape, as long as we split 
the D^mii" 1 term. 

' -A A mu m + V A m P m = F + V A mG in n 
V Am -u m = G in Q 

( p rnj _ B^ m ti m )AA m = H on S 

u m = on E,, 



(2.3.18) 



In the following, we will prove an improved estimate for (u m ,p m ) in terms of ||n m ||^fe+i/2 ( - s - ) . 
We divide the proof into several steps. 
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Step 1: Preliminaries 

To abuse the notation, within this bootstrapping procedure, we always use (u,p,r)) instead 
of (u m ,p m ,T] m ) to make the expression succinct, but in fact they should be understood as 
the approximate sequence. Also it is easy to see the term V^G will not affect the shape of 
estimate because || V A mG\\ H k-i < ||?? ,n ||^fc+i/2( S ) ||G||^fc. Hence, we still write F here to indicate 
the forcing term in first equation. 

We write explicitly each terms in above equations, which will be hired in the following. 
(2.3.19) 

dum + d 22 ui + (1 + A 2 + B 2 )K 2 d 33 ui - 2AKd n m - 2BKd 23 m 

+ {AKd 3 {AK) + BKd 3 (BK) - di(AK) - d 2 (BK) + Kd 3 K)d 3Ul + d lP - AKd 3 p = F 1 

(2.3.20) 

dnu 2 + d 22 u 2 + (1 + A 2 + B 2 )K 2 d 33 u 2 - 2AKd l3 u 2 - 2BKd 23 u 2 

+(AKd 3 {AK) + BKd 3 (BK) - di{AK) - d 2 (BK) + Kd 3 K)d 3 u 2 + d 2 p - BKd 3 p = F 2 

(2.3.21) 

dnu 3 + d 22 u 3 + (1 + A 2 + B 2 )K 2 d 33 u 3 - 2AKd l3 u 3 - 2BKd 23 u 3 

+{AKd 3 {AK) + BKd 3 (BK) - di(AK) - d 2 (BK) + Kd 3 K)d 3 u 3 + Kd 3 p = F 3 

(2.3.22) dim - AKd 3 m + d 2 u 2 - BKd 3 u 2 + Kd 3 u 3 = G 

where for all above, A, B and K should be understood in terms of n m . For convenience, we 
define 

(2.3.23) Z = C( m )P(r,) ( H^H^-i + \\G\\ 2 Hk + ||iJ||^ fc _ 1/2(s) \ 

whereC(r^o) isaconstant depending on ||%|lijfc+i/2(s) and P(??) is the polynomial of ||r/||^ fe+ i/ 2 ( S ). 



Step 2: r = k case 

For k — 1 order elliptic estimate, we have 



(2.3.24) 



M|^fc-l + ||p||#k-2 < C(l]o) ^ ||F||^- fc _3 + ||G||#fc-2 + ||iJ||^-fc_ B /2^ ^ < Z 



By lemma 2.17, the bounding constant C(r/o) only depend on ||7/o||#k+i/2. 
For i = 1,2, since (d{U,dip) satisfies the equation 



(2.3.25) 



where 



-A^u) + V A (d iP ) = F in n 

V A • (diii) = G in ft 

{{d iP )I - A (d iU ))M = H on T 

diU = on Tii 



F = 8iF + V diA ■ V A u + V A ■ V diA u - V diA p 
G = d t G- V 9iA ■ u 

H = diH - (pi - B A u)diM + O diA uAf 



Employing k — 1 order elliptic estimate, we have 



(2.3.26) 



I dn 



-i + \\d iP \\ 2 Hk -2 < Cfa>) \\F 



\jjk-3 + ||G||^fc_ 2 + \\H \\ fjk-5/2(J2) 



Since except for the derivatives of F, G and H, all the other terms on the right hand side has 
the form \\A ■ B\\ Hr , in which A = d a fj and B = d^u or d^p. These kinds of estimates can be 
achieved by lemma A.l. Because this lemma will be repeated used in the following estimates, 
we will not mention it every time and all of the following estimate can be derived in the same 
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fashion. Hence, we have the forcing estimate 



l^lli?*- 3 ll^ll H k ~ 2 ll^ll_ff fe - 5 / 2 



(2.3.27) 



< \\F\ 



+ IIGI 



-i + \\H\ 



H k- 



3/ 2(s) + C(r]o)P{v) [ \\u 



\ H k-i 



+ 



< z 

In detail, this means 

(2.3.28) \\diuf H k-i + \\d lP f H k- 2 + ||^«||h*-i + II^IIh*-* £ z 

Hence, most parts in + has been covered by this estimate, except those with 

highest order derivative of $3. 

Multiplying A to (2.3.21) and adding it to (2.3.19) will eliminate the d 3 p term and get 

{d llUl + Ad u u 3 ) + {d 22 ui + Ad 22 u 3 ) + {I + A 2 + B 2 )K 2 {d 33Ul + Ad 33 u 3 ) - 2AK(d 13Ul + Ad 13 u 3 ) 
-2BK(d 23Ul + Ad 23 u 3 ) + (AKd 3 {AK) + BKd 3 (BK) - d^AK) - 8 2 {BK) + Kd 3 K){d lUl + Ad lU3 ) 
+dip = F 1 + AF 3 

Then taking derivative d k ~ 2 on both sides and focus in the term (l+A 2 + B 2 )K 2 (d^ui+Ad^u 3 ), 
the estimate of all the other terms in H° norm implies that 



(2.3.29) 

Similarly, we have 
(2.3.30) 



d 3 u\ + Ad 3 u 3 



ho 



< Z 



< Z 



dlu 2 + Bdlu 3 

Rearrange the terms in (2.3.22), we get 

(2.3.31) K{1 + A 2 + B 2 )d 3 u 3 = G — d lUl - d 2 u 2 + AK(d 3Ul + Ad 3 u 3 ) + BK{d 3 u 2 + Bd 3 u 3 ) 

Taking derivative c^ -1 on both sides, focusing in the term K(l + A 2 + B 2 )d 3 u 3 employing all 
of the estimate we have known, we can show 



(2.3.32) 



d k 3 u 3 



H° 



<Z 



Combining (2.3.29), (2.3.30) and (2.3.32), it is easy to see we can get 



(2.3.33) 



H° 



+ 



d k 3 u 2 



HO 



< Z 



Plugging this to (2.3.21) and taking derivative <9 3 on both sides, we get 
(2-3.34) \\ d "~ lp \\ 2 H0 



< Z 



Combining this with all above estimate, we have proved 



(2.3.35) \\uf H k + 

Therefore, we have proved the case r = k. 



\H k 



1 <z 



Step 3: r = k + 1 case: first loop 

For i,j = 1,2, since (dijU,dijp) satisfies the equation 



(2.3.36) 



-A A (diju) + V A (dijp) = F in tt 

V 'a " (diju) = G in CI 

((d ijP )I-B A (d ijU ))M = H on S 

, (diju) = on S 6 
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where 

F = dijF + V dijA ■ V A u + X7 A • V dijA u + V diA ■ V djA u + V 9jA ■ V 9i AU + V diA ■ V A {d 3 u) 

+VdjA ■ ^A(diu) + V A ■ V diA (dju) + V A ■ V d]A {diu) - V 9ijA p - ^d,A(djP) - ^d,A(diP) 
G = dijG -V 9ijA -u-V diA - (dju) -V 9jA - (diu) 

h = d l3 H - {pi - o^dijM - {{d jP )i - B A (d jU ) - B 9jA u)d t M - {{d iP )i - o A (d t u) - n diA u)d 3 N 

+ (O dtjA u + B diA (d jU ) + B djA (d iU ))M 
Employing k — 1 order elliptic estimate, we have 

(2.3.37) ll&iiillifc-i + Wdnpf^-i < F 2 + G 2 + H 2 <Z 



H<=-5/2(E) 



where the forcing estimate can be taken in a similar argument as in r = k case. 

In detail, this is actually 

(2.3.38) 

||diiu||^ fc -i + \\di 2 u\\ 2 Hk -i + \\d 2 2u\\ 2 H k-i + \\dnp\\ 2 H k-2 + \\d 12 p\\ 2 H k-2 + 11^22^11^-2 ~ z 



Step 4: r = k + 1 case: second loop 

Similar to r = k case argument, for i = 1,2, taking derivative d^~ 2 di on both sides of (2.3.29) 
and focus in the term <9f <9jiii + Ad^diU 3 , we get 



(2.3.39) 

Similarly, we have 
(2.3.40) 



dld iUl + Adld iU3 



did iU2 + Bd$d iU3 



HO 



<z 



H° 



<z 



Plugging in this result to (2.3.22) and taking derivative d% 2 di on both sides, we will get 



(2.3.41) 



dld lU3 



HO 



<z 



An easy estimate for these three terms implies 



(2.3.42) 



d$di 



ii 



HO 



< z 



Plugging this to (2.3.21) and taking derivative <9g 2 di on both sides, we get 



(2.3.43) 



it^P 



HO 



<z 



Combining all above estimate with estimate of previous steps, we have shown 

(2.3.44) \\d3iu\\Hk-i + \\d*P\\%k-* 
In detail, this is actually 

(2.3.45) Insulin*-! + \\d2M\l^ + WdiwWl*-* + II^pIIh*-" £ z 



Step 5: r = k + 1 case: third loop 

Again we use the same trick as above. Taking derivative d^ 1 on both sides of (2.3.29) and 



focusing in the term dn +1 u\ + Adn +1 u 3 , we can bound dn +1 u\ 



Bd^ +1 u 3 . 



A8l +1 u ?! . 



Similarly, we control 



Plugging in this result to (2.3.22) and taking derivative 1 on both sides, 



we can estimate dl +1 u 3 . Then we have 
(2.3.46) 



HO 



< Z 
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Plugging this to (2.3.21) and taking derivative d 3 on both sides, we get 
(2.3.47) 



d k p 



H° 



< z 



Combining this will all above estimate, we get 
(2.3.48) 



||^3«Hh*-i + \\d33P\\l k -2 < Z 



Step 6: r = k + 1 case: conclusion 

To synthesize, (2.3.38), (2.3.45) and (2.3.48) imply that all the second order derivative of u in 
k — 1 norm and p in k — 2 norm is controlled, so we naturally have the estimate 

(2.3.49) + IIpIIh* £ c (vo)P(v) ( ll^ll^-i + \\G\\ 2 h h + ||^||^-i/ 2 ) 

This is just what we desired. 

Now let us go back to original notation, since P(-) is a fixed polynomial, this gives an estimate 
\\u m \\W + \\p m \\ 2 H* < C{ m )P{n rn ) ( + \\G\\ 2 Hk + ||F||^ fc _ 1/2(s) ) 

<C( Vo )P(v)( \\P\\ 2 Hk -i + \\G\\ 2 Hk + \\H\\ 

(2.3.50) 



<C( Vo )P( m )[ \\F\\ 2 Hk - 1 + \\G\\ 2 Hk + \\H" 2 



Hfc-i/2(E) 



< C{vo) ^ ll-^lllffc-i + \\G\\ 2 H k + \\H\\ 2 H k-i/2^) ^ 



where C(rj ) depends on WvoWh^ 1 / 2 ^)- 

This bound implies that the sequence (u m ,p m ) is uniformly bounded in H k+1 x H k , so we can 
extract weakly convergent subsequence u m u° and p m p°. 

In the second equation of (2.3.16), we multiplied both sides by J m w for w £ to see that 

f GwJ m = f (V Am ■ u m )wJ m = - [ u m - {V A m W )J m -> - / u° ■ (V A w)J = f (V A ■ u°)wJ 
Jn Jn Jn Jn Jn 

(2.3.51) 

which implies ■ u° = G. Then we multiply the first equation by wJ m for w £ W and 
integrate by parts to see that 

(2.3.52) I [ A m U m :B A m W J m - [ p m V Am -w,J m = [ F ■ w,J m — f H w 
2 Jn Jn Jn Jv. 

Passing to the limit m — > oo, we deduce that 

(2.3.53) - / B A u° : B A wJ - [ p°V A ■ wJ = [ F ■ wJ — [ H-w 
2 Jn Jn Jn is 

which reveals, upon integrating by parts again, (u°,p°) satisfies (2.3.1). Since (u,p) is the unique 
strong solution to (2.3.1), we have u° = u and p° = p. This, weakly lower semi-continuity and 
estimate (2.3.50) imply (2.3.54). □ 

Remark 2.19. The key part of this proposition is that as long as we can deduce r/o £ H k+1 / 2 (Y>) 
and 7] £ H k+l / 2 {T,), we can achieve u £ H k+1 and p £ H k , which is the highest possible 
regularity we can expect. 

Remark 2.20. By our notation rule, since C(r/o) depends on and is given implicitly, we 
can take it as a universal constant, so the estimate may be rewritten as follows. 

(2.3.54) IMIfp + IIpH^-i ^5 ll-^1l_f/ r - 2 + ll^ll/f-i + 11-^1 3/ 2 (s) 
for r = 2, . . . , k + 1 
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2.3.2. A-Poisson Equation. We consider the elliptic problem 

= f in Q 
p = g on £ 

VaP • v = h on E5 

For the weak formulation, we suppose / G V*, g G # 1/2 (£) and ft G Jf-V2(£ 6 ). Let p G -ff 1 
be an extension of 5 such that supp(p) C {— inf(b)/2 < X3 < 0}. Then we can switch the 
unknowns to q = p — p. Hence, the weak formulation of (2.3.55) is 

(2.3.56) {V A q, V A <P)n° = ~(^AP, V A <t>) u o - (/, <f>)v + {h, (f>)- 1/2 

for all <fi G V, where {-,-)v* denotes the dual pairing with V and (■,■)— 1/2 denotes the dual 
pairing with iJ -1 / 2 (£;,). The existence and uniqueness of a solution to (2.3.56) is given by 
standard argument for elliptic equation. 

If / has a more specific fashion, we can rewrite the equation to accommodate the structure of 
/. Suppose the action of / on an element 4> G V is given by 

(2.3.57) (/, 4>) v * = (fo, <P)u° + (Fo, Va<P)u° 
where (f ,F ) G H° x H° with \\f Q \\ 2 H0 + || 

f o\\ho = WfWv*- Then we rewrite (2.3.56) into 

(2.3.58) (V A P + ^0, V40) w o = -{fo, 4>) w + (h, <f>)_ 1/2 
Hence, it is possible to say p is a weak solution to equation 

f V4 • (VaP + F ) = fo 

(2.3.59) I P = 9 

{ (V A P + F )-u = h 

This formulation will be used to construct the higher order initial conditions in later sections. 

Lemma 2.21. Suppose that 77 G H k+1 / 2 (Y>) for k > 3 such that the map $ is a C 1 diffeomor- 
phism ofVL to Jl' = $(fj). /// G H°,g G H 3 / 2 and h G H 1 / 2 , then the equation 2.3.55 admits a 
unique strong solution p G H 2 . Moreover, for r = 2, . . . , k — 1, we have the estimate 

(2.3.60) \\p\\ Hr < C(n) ( \\f\\ Hr -2 + ||5ll^-i/2 (s) + ll^ll^-3/ 2(S6) \ 

whenever the right hand side is finite, where C{n) is a constant depending on \\r]\\ H k+i/2^)- 

Proof. This is exactly the same as lemma 3.8 in [1], so we omit the proof here. □ 

Next, we will prove the bounding constant for the estimate can actually only depend on the 
initial surface. Since we do not need optimal regularity result for this equation, we do not need 
the bootstrapping argument now. 

Proposition 2.22. Let k > 6 be an integer suppose that r\ G H k+1 / 2 (Y>) and r/o G H k+1 / 2 (Y,). 
There exists eo > such that if \\rj — ??o || ij fc-3/2 < eo, then solution to 2.3.55 satisfies 

(2-3.61) ||p|| Hr < C(r]o) ( \\f\\ H r-2 + 115-11^-1/2(2) + I \ h \ lff-3/2(s 6 ) 1 

for r = 2, . . . , k — 1, whenever the right hand side is finite, where C(??o) is a constant depending 
on Iholl/ffc+i/a^)- 

Proof. The proof is similar to lemma 2.17. We rewrite the problem as a perturbation of the 
initial status. 

' A Ai) p m = f + f m in ft 
(2.3.62) < p m = g + g m on E 

, ^A () P m -v = h + h m on S 6 

The constant in this elliptic estimate only depends on the initial free surface. We may estimate 
f m , g m and h m in terms of p m , r/o and £ m = n m — r/o- The smallness of £ m implies that we can 
absorb these terms into left hand side, which is actually (2.3.63). □ 
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Remark 2.23. Similarly, we can also take this constant C(j)q) as a universal constant and 
rewrite the estimate as follows. 

(2.3.63) \\p\\H r ~ II/IIh'- 2 + II5'IIh'- 1 /2(e) + IHIiF-3/2(s 6 ) 

for r = 2, . . . , k — 1. 

2.4. Linear Estimates. 



(2.4.1) 



dtu — A A u + V ' A u = F in £1 

V A ■ u = in Q 

(pi - O a u)M = H on £ 

u = on £5 

In this section, we will study the linear Navier-Stokes problem (2.4.1). Following the path of 
I. Tice and Y. Guo in [1]. We will employ two notions of solution: weak and strong. Then we 
prove the wellposedness theorem from lower regularity to higher regularity. To note that, most 
parts of the results and proofs here are identical to section 4 of [1], so we won's give all the 
details. The only difference is the bounding constants related to free surface since our surface 
has less restrictions. 

In this section, we always assume P(z) as a polynomial of z. 

2.4.1. Weak Solution and Strong Solution. The weak formulation of linear problem (2.4.1) is as 
follows. 

(2.4.2) (d t u,v) C 2 H o + - (B A u,B A v) C 2 H o - (p,V A ■ v) C 2 H o = (F,v) C 2 H o - (H,v) C 2 H o {S) 

Definition 2.24. Suppose that u G H°, F G W* and H G L 2 ( [0, T] ; H~ x / 2 (£)) . If there exists 
a pair (u,p) achieving the initial data uq and satisfies u G W, p G L 2 ([0, T]; H°) and dtu G W* ; 
such that the (2.4-2) holds for any v G W, we call it a weak solution. 

If further restricting the test function v G X, we have a pressureless weak formulation. 

(2.4.3) {d t u,v) C 2 H o + ^(\D A u,B A v) C 2 n o = (F,v) C 2 H o - (H,v) C 2 n o (s) 

Definition 2.25. Suppose that u G y(0), F G X* and H G L 2 ([0, T]; iJ _1 / 2 (S)). If there 
exists a function u achieving the initial data u$ and satisfies u G X and dtu G X* , such that 
the (2.4.3) holds for any v G X, we call it a pressureless weak solution. 

Remark 2.26. It is noticeable that in O, W* C X* and ||u||^» < ||u||yy»; on S ; for u G X , we 
have u\y, G L 2 ([0,T]; H 1 ^ 2 (T,)) . 

Since our main aim is to prove higher regularity of the linear problem, so the weak solution 
is a natural byproduct of the proof for strong solution in next subsection. Hence, we will not 
give the existence proof here and only present the uniqueness. 

Lemma 2.27. Suppose u is a pressureless weak solution of (2.4-1) ■ Then for a.e. t G [0, T], 
(2.4.4) 



\ [ JHt)\ 2 +\ f [ J|P^| 2 = ^ / J(0)N| 2 +i f [ \u\ 2 d t J+f [ J{F.u)-f [ H- 
z Jn z Jo Jn z Jn z Jo Jn Jo Jn Jo is 

Also 

(2.4.5) \\u\\ 2 L oo H o + \\u\\ 2 L 2 H i < C {v) ^ ||«o||^o + \\F\\x* + |lia H -i/2 (E) ^ 



u 



where 



(2.4.6) C (n) = P(l + \\n\\ H5/ 2 {S) + ||^|| H5/2(S) ) exp (tP(1 + \\d tV \\ H5/ 2 { ^ 

Proof. The proof is almost the same as that of lemma 4.1 in [1]. The only difference is that we 
should replace the lemma 2.1 in [1] used in that proof by our lemma 2.5 now. □ 
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Proposition 2.28. // a pressureless weak solution exists, then it is unique. 

Proof. If u 1 and u 2 are both pressureless weak solution to (2.4.1), then w = u 1 — u 2 is a 
pressureless weak solution with F = H = w(0) = 0. Then the estimate (2.4.5) implies that 
w = 0. Hence the solution is unique. □ 

Next, we define the strong solution. 
Definition 2.29. Suppose that 
uo G H 2 n X(0) 

(2.4.7) F G L 2 ([0, T]; i? 1 ) D C°([0, T];H°) d t F G X* 

H G L 2 ([0,T};H 3 / 2 (Z))nC°([0,T];H l / 2 (Z)) 8 t H G L 2 ([0, T]; H' 1 / 2 ^)) 

If there exists a pair (u, p) achieving the initial data uq and satisfies 

uGL 2 ([0,T];tf 3 )nC< ([0,T];tf 2 )n* d t u G L 2 ([0, T}; tf 1 ) n C°([0, T]; H°) 
{Z ^> 8 2 u£X* pG L 2 ([0,T];^ 2 )nC°([0,T];^ 1 ) 

smc/i £/ia£ i/iey satisfies (2.4-1) in the strong sense, we call it a strong solution. 

2.4.2. Lower Regularity Theorem. 

Theorem 2.30. Assume the initial data and forcing terms satisfies the condition (2.4-7). De- 
fine the divergence-A preserving operator D t u by 

(2.4.9) D t u = d t u-Ru for R = d t MM~ l 

where M is defined in (2.2.30). Furthermore, define the orthogonal projection onto the tangent 
space of the surface {X3 = tjq} according to 

(2.4.10) n (v) = v- (vAf )Af \Af \- 2 for Af = (-dir) , -d 2 r) , 1) 
Suppose the initial data satisfy the compatible condition 

(2.4.11) IL (H (0) + D Ao u M ) = 

Then the equation (2.4-1) admits a unique strong solution (u,p) which satisfies the estimate 
(2.4.12) 

+ \\d 2 u\\ 2 w * + |M| 

< L(rj) ( Huoll^ + 11^(0)11^0 + ||^(0)||^i /a + \\F\\\ 2m + \\d t F\\ 2 x * + \\H\\ 2 l2hZ/2 + \\d t H\\ 2 L2H _ 1/2 
where 

(2.4.13) L(r?) = P(l + K{rj)) exp (tP{1 + K{ V ))^j 
for 

(2.4.14) K(rj) = suV t ( ||r?(t)||^9/ 2(s) + \\d t v(t)\\ H V H i:) + 

The initial pressure p(Q) G H 1 , is determined in terms of u$, t]q, F(0) and H(0) as the weak 
solution to 

( V Ao -(V Ao p(0)-F(0)) = -V Ao (R(0)uo) in n 

(2.4.15) I p(0) = (H(0)+B Ao u N )-N \No\- 2 on £ 
( (V Ao p(0)-F(0))-v = A Ao u -v on S 6 

in the sense of (2.3.59). 

Also, D t u(0) = dtu(0) — R(0)uo satisfies 

(2.4.16) D t u(0) = A Ao u - V Ao p(0) + F(0) - R(0)u o G y(0). 
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(2.4.17) 



Moreover, (D t u, dtp) satisfies 

' d t {D t u) - A A (D t u) + V A (d tP ) = D t F + G F in fi 

V.4 • (An) = in tt 

S A (d tP ,D t u)N = d t H + G H on S 

D t u = on S;, 

in the weak sense of (2.4-3), where 

G F = -(R + d t JK)A A u-d t Ru+(dtJK + R + R T )V A p + V A -(B A (Ru) 
G H = B A (Ru)M - (pi -B A u)d t M + B dtA uM 



RB A u + n dtA u) 



Proof. Since this theorem is almost identical to theorem 4.3 in [1], we will only point out the 
differences here without giving details. In order for contrast, we use the same notation and 
statement as [1] here. 

(1) The only difference in the assumption part is that in theorem 4.3 of [1], K(rj) is suf- 
ficiently small, but our K(n) can be arbitrary. Hence, in [1], we can always bound 
P(l + K(n)) < 1 + K(n), but now we have to keep this term as a polynomial (this is 
the only difference of L(rj) between these two theorems). 

(2) It is noticeable that each lemma used in proving theorem 4.3 of [1] has been recovered 
or slightly modified in our preliminary section. To note that these kind of modification 
merely involves the bounding constant due to the arbitrary K(rj), but does not change 
the conclusion and the shape of estimate at all. Therefore, all these modification will 
only contribute to P(l + K(rj)) term. 

(3) In elliptic estimate, we now use proposition 2.18 to replace proposition 3.7 in [1], which 
gives exactly the same estimates. 

Based on all above, we can easily repeat the proving process in [1] with no more work, so our 
result naturally follows. □ 

2.4.3. Initial Data and Compatible Condition. Before starting to discuss the higher regularity 
result, we first need to determine the initial data to higher order derivatives and the correspond- 
ing compatible conditions. Although this part is also identical to section 4.3 of [1], for clarity, 
we repeat this process here. 
Define a vector field on Q: 

(2.4.18) g°{F,v,q) = A A v-V A q + F-Rv 
Define function f on fl, g on S and f) on S5 according to 

f(F,v) = V A -(F-Rv) 

(2.4.19) g(H, v) = (H + B A vAf) ■ N \N\~ 2 

t,(F,v) = (F + A A v) ■ v 

Define the quantities related to the initial data: 

N n N-l 



(2.4.20) 
(2.4.21) 



Ho(u,p) 
fC (u ) 



j=0 

11 1 1 2 
\\ u o\\h 2N 



i=o 



(2.4.22) 
(2.4.23) 

(2.4.24) 



n (v) 

Ko{F,H) 



N 

ii?ko)ii;U + i/ 2 (£) + E||^ ) 

11 n 2 

||%||#2iV+l/2( S ) 

N-l 2 N-l 

y |k'F(o) + y \\d 3 t H(o) 

3=0 3=0 



^2JV-2j+3/2( S ) 



#2JV-2j-3/2(£) 
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Furthermore, we need to define mappings Q\ on Q and Q2 on S. 
(2.4.25) 

g 1 ^, g) = -(i? + 3 t JK)A^ - ftifo, + (d t + R + i2 T )V^g + ■ (p A (Rv) - RO A v + B dtA v) 
G 2 (v, q) = O a (Rv)M - (ql - A v)d t M + B dtA vM 

The mappings above allow us to define the arbitrary order forcing terms. For j = 1, . . . , N, we 
have 

F° = F 
H° = H 

(2.4.26) Fj = DtF j-i + = ^ F + Dlg^Di-^u, dl~ l ' l p) 

H i = DtH i-i + g\ D {-\ di-'p) = d{h + n=o ^(iT' - V ar' _1 P) 

Finally, we define the general quantities we need to estimate as follows. 

N , N-l 

(2.4.27) = 

3=0 
N 

(2.4.28) 2?(«,p) = 

3=0 

(2.4.29) /C(u,p) = S(u,p) +V(u,p) 



X jX ff2N-2j 



+ E Wp 



3=0 



£ J ooff2N^2j-l 

N-l 



L 2 H 2N-2j + l 



+ 



+ E H p 

3=0 



L 2 H 2N-2j 



N 



(2.4.30) £(77) 

(2.4.31) V{r]) 



E in 1 

3=1 



L oo#2iV-2j+3/2( S ) 



JV+1 



||r/||^2 H 2iv+i/2( E ) + \\d t r]\\ 2 L 2 H 2N-i/2^ + E 

3=2 

(2.4.32) /Cfa) = £(»/)+ 2>fa) 



L 2#2JV-2j+5/2( S ) 



(2.4.33) K{F,H) 



N-l 

E ||3> 

AT-1 

+ E|h F 



N 



L 2 H 2N-2j 



, + ll«,+Ep 

AT-1 



L 2^2JV-2j-l/2( S ) 



Looff2N-2j-2 



+ E N F 



J=0 



L oo/f2iV-2j-3/2( S ) 



In the following, we will present some preliminary results to estimate the forcing terms and 
initial data, which will be used both in constructing the higher order initial data and proving 
the higher regularity theorem. 

Lemma 2.31. For m = 1, . . . , N — 1 and j = 1, . . . , m, the following estimates hold whenever 
the right-hand sides are finite: 

(2.4.34) 



\ FJ \ 

3-1 

1=0 



L 2 H 2m,-2j + 



1 II mil < 

1 ^ 11-" llL 2 H 2m - 2 j+3/ 2 (S) ~ 



P(1 + )C( V ))()C(F,H) 





2 




2 


\d\u 


+ 








£2^2m-2j+3 




]^oofj2m-2j+2 



3-1 



+J2\\ d tp 



1=0 



l2u2m-2j+2 



+ 



d l tP 



J^ao I]2m-2j + l 



(2.4.35) 



I V3 1 1 ^ 4-11 fT 3 11^ < 

1^ ||^oo^2m-2i ' W 12 ||^oo^2m-2i+l/2( S ) ^ 



( II 

^ 1=0 



a 



Lcx>fl2m,-2j + 2 



+ 



d l tP 



£00 Jj2m-2j+l 
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\d t F m \\ x * + \\dtH m \\ 2 L2H - 1/2 ^ 



< P(l + /C(r?)) (/C(P, P) + \\dTu\\\ 2H2 + ||5 t >||i 2f , 1 



m— 1 

+ TM 



1=0 



u 



L^H 2 



+ 



m—1 



L 2 H 3 Z—/ II 
Z=0 



+ 



d l tP 



L 2 H 2 



Proof. The estimates follow from simple but lengthy computation, involving standard argument. 
Since this lemma is almost the same as lemma 4.5 in [1] and the argument are exactly the same, 
we will only present a sketch of the proof here. 

After we expand all the terms with Leibniz rule and plug in the definition of F 3 and H 3 , it is easy 

1 1 2 1 1 2 

to see that the basic form of estimate is ||Xy|| i2 ^fc or || AY))^^ where X includes the terms 
only involving f) and Y includes the terms involving u, p, F or H. Employing lemma A.l, we can 

2 2 2 2 2 2 

always bound \\XY\\ L 2 Hk < \\X\\ L2[{kl \\Y\\ L<xHk2 or \\XY\\ L 2 Hk < ||X|| L<x>Hfcl \\Y\\ L2Hk2 , and 
also HXFll^oo^fc < H-X'Hiootffci H^ll^so^ • The principle of choice in these two forms for L 2 H h 
is not to exceed the order of derivative in the definition of right-hand sides in the estimates. A 
detailed estimate can show this goal can always be achieved. Hence, we can show (2.4.34) and 
(2.4.35). 

For the estimate (2.4.36), we always have the trivial bound 
(2-4.37) \\f\\%. < \\ff L2H0 

In dtF m , except the term d^F, we can always use the same argument as above, so our estimate 
can be easily shown. □ 

Lemma 2.32. For j = 1, . . . , N - 1, 

(2.4.38) 

<P(l+H (r]))(lC (F,H) 



j'-i 



1=0 



+ 



d l tP (0) 



]-[2N-2l-l 



Proof. It is similar to proof of lemma 2.31, so we omit it here. 
Lemma 2.33. If k = 0, . . . , 2N — 1 and v is sufficiently regular, then 
(2-4.39) \\d t v - D t v\\ 2 L 2 Hk < P(l + K{n)) \\v\\ 2 L 2 Hk 

and if k = 0, . . . , 2N — 2 and v is sufficiently regular, then 
(2-4.40) \\d t v - D t v\\ 2 LooHk < P(l + K{n)) \\v\\ 2 LooHk 

, N — 1 and j = 1, . . . , m, assuming v is sufficiently regular, then 

<P(l + JC( V )) 



Ifm = l, 
(2.4.41) 



cP t v — D\v 



l2u2m-2j+Z 
3-1 

E 

1=0 



\dlv 



l2fj2m-2j+3 



+ 



l i <xfj2m-2j+2 



□ 



(2.4.42) 
and 

(2.4.43) 



d 3 t v 



X J ooJj2m-2j + 2 



i-i 

<p(i+ic( V ))J2 

1=0 



L<x f]2m-2j+2 



\\d? +1 v - d t D?v\\ 2 L2Hl + \\dT +2 v - d 2 t D?v\\ 2 Xf < P(l + Kb)) 



EN 



1=0 



L 2 m 



+ 



d l t u 



L°°H 2 



+ m 



m+l 



X* 
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Proof. In a similar fashion as proof of lemma 2.31, we can easily show this. □ 
Lemma 2.34. If j = 0, . . . , N and v is sufficiently regular, then 



(2.4.44) div(0) - Div(0) * < P(l + H (r,)) £ \\d l t v(0) 



H 2N- 

1=0 



2 



JJ2N-21 

Proof. The proof is similar to lemma 2.31, so we omit it here. □ 



Lemma 2.35. Suppose that v, q, G\ and G2 are evaluated at t = and sufficiently regular for 
the right-hand sides of the following estimates to make sense. For j = 0, . . . , N — I, we have 

(2.4.45) \\G°(G 1 ,v,q)\\ 2 H2N _ 2j _ 2 < P(l + ||»/(0)||^ (E) + ||^(0)||^_i (E) ) 

^ ||f||^2iV-2j + ||oJ| jj2N-2j-l + || G\ || 2j— 2 

Ifj = 0,...,N-2, then 

(2.4.46) ||f(Gi,«)||^ aJV -«-3 + ||0(G 2 , V )||^_ 2 ,_3 /2(E) + ||f)(Gi , v)\\ 2 H 2N-2j-5/2( Ti ) 

< P(l + ||?7(0)||#2JV( 2 )) ^ ||Gl||^-2JV-2j-2 + ||G2||^2iV-2j-3/2( S ) + || 1> || #2JV-2j ^ 

For j = N — 1, if V^ • t) = in fi, t/ien 

(2.4.47) ||fl(G2,«)||HV»(E) + llb(Gi,«)|| <P(1 + N(0)|| 

|Gl||^2 + ||G 2 ||^1/2( S ) + |M|#2 

Proof. The proof follows from a similar fashion as those of lemma 2.31. The only thing to note 
is that V^ ■ v = implies V^ () • A^ Q v = 0, so we have \\Aa v ' ^llff-V2(E 6 ) ~ II^Aj^II^o- 1=1 

Now we start to determine the initial data. We assume that uq G H 2N (Q), 7/0 G # 2iV+1 / 2 (£) 
and JCq(F, H) < 00. In the following, we will iteratively construct D J t u(0) for j = 0, . . . , N and 
5£p(0) for j = 0,...,N- 1. 

For j = 0, F°(0) = F(0) G H 2N ~ 2 , H°(Q) = H(0) G P 2iV - 3 / 2 (£) and £>°u = u G F 2iV . 
Suppose now F'(0) G H 2 "- 21 ' 2 , H l (Q) G #2Af-2Z-3/2 (S ) and d j u ( ) £ R 2N-2l are given 

for < I < j G [0,N - 2], we will define 5£p(0) G H 2N - 2 i~ l as well as r^' +1 u G H 2N - 2 \ 
Fi +1 (0) G H 2N - 2 i~ 4 and i^' +1 (0) G H 2N -^- 7 / 2 {T,). By virtue of (2.4.46), we know that / = 
f{Fi(0),Diu(0)) G F 27V - 2 J- 3 , 5 = 0(^(O),^u(O)) G H 2N ' 2 ^' 3 / 2 {^) and /i = t)(F^(0), D j t u{Q)) G 
P 2 ^ -2 -? -5 / 2 (£;,). This allows us to define dfp(0) as the solution to (2.3.55) with this choice of 
/, g, h and then the proposition 2.22 with k = 2N implies df.p(0) G iJ 2iV_2 J _1 . The estimate 
(2.4.38), (2.4.44) and (2.4.45) allow us to define 

D{ +1 u{0) = g°{Fi(0),Diu{0),d{p{0)) G H 2N - 2 ^ 2 

(2.4.48) FJ +1 (0) = AP j (0) + g\D{u(ti), %p(0)) G H 2N - 2 ^ 4 
Hi +1 (0) = d t W{0)+g 2 {Diu(0),d{p{0)) G i? 2JV - 2 ^ 7 / 2 (S) 

Iteratively, we can construct all the data except D^u(0) and dj*~ 1 p(Q). 

In order to complete this construction, we must enforce our compatible condition for j = 
0, . . . , N — 1. We say that the compatible condition is satisfied if 



(2.4.49) 



D}u(0) G H 2 n X(0) 
n o (^(0)+B A)J D^(0)AA ) = 



It is easy to see our construction satisfies this compatible condition for all j, so V_4 () -Df 1 u(0) = 
0, which means we can use estimate (2.4.47) to see that g = g(F Ar " 1 (0), _1 n(0)) G # 1/2 (£) 
and h = f}(H N - 1 (0),Df- 1 u(Q)) G H' 1 / 2 ^). We can also define /„ = -V^ o (i?(0)Df- 1 «(0)) 
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with /o G H° and F = -F N - 1 (0) G Hence, we can define <9f~V ) e i^ 1 as the weak 
solution of (2.3.55) in the sense of (2.3.59). Then we can just define 

(2.4.50) D?u{0) = g°(F N - 1 (0),D^ 1 u{0), 8^-^(0)) G H° 

In fact, above construction guarantees that D^u(0) G 3^(0). Furthermore, the bounds (2.4.38), 
(2.4.45) and (2.4.46) provides us with the estimate 

(2.4.51) 

N II 2 N ~ 1 ll 2 / \ 

T.\\ D lM H ^ + ^ ¥ m mx-v-i Z P ( 1 + n ^))( \\uof H , N +)C (F,H)) 

3=0 j=0 V 7 

Owing to estimate (2.4.44), the above estimate also holds when replace D\u with &[u, i.e. 

(2.4.52) n (u, P ) < P(l + H (v)) ( ho\\ 2 H 2N + JC (F, H^J 
2.4.4. Higher Regularity Theorem. 

Theorem 2.36. Suppose the initial data uq G H 2N , t]q G H 2N+1 / 2 (Y,) and the initial condition 
is constructed as above to satisfy the j th compatible condition for j = 0, 1, . . . , N — I, with the 
forcing function K.(F,H) + K,q{F,H) < oo. Then there exists a universal constant To(n) > 
depending on C(n) as defined in the following, such that if < T < To(n), then there exists a 
unique strong solution (u,p) to the equation (2.4-1) on [0, T] such that 

d{u G L 2 ([0, T]; H 2N ~ 2j+l ) n C°([0, T]; H 2N ~ 2j ) for j = 0,...,N 

(2.4.53) d> t p£ L 2 {[{),T]-H 2N - 2 ^r\C Q {^,T] ] H 2N - 2 3- 1 ) for j = 0,...,N-l 
d? +1 u G X* 

Also the pair (u,p) satisfies the estimate 

(2.4.54) K{u,p) <£( V )^\\u \\ 2 H2N +}C (F,H) + ]C(F,H)^ 
where 

(2.4.55) C(n) = P(l + U {n) + K{rj)) exp ( TP{1 + K[rf)) 



Furthermore, the pair (D^u^d^p) satisfies the equation 

' d t (p\u) - A A (D J t u) + V a (Hp) = FJ in n 



(2.4.56) 



V A - (D J t u) = in n 



S A {d> t p,D{u)N = W on E 

D j t u = on E fe 

in the strong sense for j = 0, . . . , N — 1 and in the weak sense for j = N . 

Proof. Similar to lower regularity case, this theorem is almost identical to theorem 4.8 in [1]. 
Due to the same reason that our free surface n can be arbitrary, the only difference here is 
that we keep the polynomial of JC(n) and Ho(n) in the final estimate which cannot be further 
simplified. Hence, our result easily follows. □ 

3. Transport Equation 
3.1. Introduction. In this section, we will prove the wellposedness of the transport equation 

dtn + u\d\n + ^2^2?? = 113 on E 



1 "(0) = % 

Throughout this section, we can always assume u is known and satisfies some bounded proper- 
ties, which will be specified in the following. 

3.2. Transport Estimates. 
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3.2.1. Initial Data and Compatible Condition. In order to derive the wellposedness theorem of 
the transport equation for arbitrary order, we should first determine the initial data of arbitrary 
order. In this section, we can always assume T-Lo(u,p) < oo and JC(r)o) < oo, then we try to 
define dfrj(0) for j = 1, . . . , N. 

First, since we know r/(0) = r/o G H 2N+1 / 2 (Y,), we can define dtr/(0) via the transport equation. 

(3.2.1) dtv(0) = u(0) ■ AT(0) G if 2JV -V2( S ) 

Notice that u(0) G H 2N ~~ 1 / 2 (E) by trace theorem and this term will dominate. In above 
estimate, we mainly employ the fact that H r (T,) is an algebra for r > 5/2. 
For j = 2, ... , N, assuming r?(0) G H 2N+1 / 2 (Y,) and d^rj(0) G H 2N - 2k+3 / 2 (2l) for k = 1, . . . ,j-l, 
then we try to determine &(r](0). Considering the transport equation dtf] = u ■ M and u is 
sufficiently regular, we can always take j — 1 temporal derivative on both sides. 

(3.2.2) d{n = £ C j ~ 1 ) d\u ■ dl- l - l M 

1=0 ^ ' 

Since d\u G C°([0, T]; H 2N ~ 2k ~ 1 / 2 (Y,)), we can evaluate the right-hand side of above equation 
in t = 0. Then 

j'-i 

(3.2.3) ^r?(0) = ^«(0) • ar 1_i AA(0) G H 2N ~ 2 ^ 2 (£) 

z=o 

This completes the construction of initial data. Furthermore, it is easy to see that we have the 
estimate 

(3-2.4) Ho(t,)<P(1 + Ho(u, P ))\\ 

3.2.2. Wellposedness of Transport Equation. Define the quantity 



v 



2 



N 2 N-l 

(3.2.5) Q{u) = Y \\&>u + Y M 

v ' v ' Z—/ II 1 jj. g2N-2j+l Z_-/ || 1 

3=0 ' 3=0 

Obviously, we have the relation Q(u) < )C(u,p). 

Theorem 3.1. Suppose that /Co(f/o) < oo and Q(u) < oo. Then the problem (3.1.1) admits a 
unique solution r] satisfying K,(q) < oo and achieving the initial data c^r/(0) for j = 0,...,N. 
Moreover, there exists a < T{u) < 1, depending on Q(u), such that ifO<T<T, then we 
have the estimate 

(3.2.6) 1C( V ) < (l + /Cofa,))P(l + Q(«)) 
Proof. We divide the proof into several steps. 

Step 1: Solvability of the equation: 

Similar to proof of theorem 5.4 in [1], based on proposition 2.1 in [3], since u G L 2 ([0, T]; H 2N+1 / 2 (T,)), 
then the equation (3.1.1) admits a unique solution n G L°°([0,T];F 27V+1 / 2 (£)) achieving the 
initial data 77(0) = 770- 

Step 2: Estimate of £ {rj): 
By lemma B.l, we have 
(3.2.7) 



||t7|| L oo #2^+1/2(2) < exp J ||w(£)||#w+i/2(e) ^3(4)11^-2^+1/ 

for C > 0. The Cauchy-Schwarz inequality implies 

(3.2.8) 11^)11^+1/2^ dt^cj \\u(t)\\ H 2 N+1 dt<Vf^/Q(u) 



1 : -p\ dt 
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So if T < 1/(2Q(«)), then 

(3.2.9) exp^C^ \\u(t)\\ H 2N+i/2 {i:) d?j < 2 

Hence, we have 

(3-2.10) h\\loo H 2N + „^ ) < /C (%) + TQ{u) < KoiVo) + 1 

Based on the equation, we further have 

(3.2.11) 115^11^00^2^-1/2(2) < ||^3|Il«>//2JV-1/2( S ) + 11^11^00^2^-1/2(2) | M | £oo^2iV+l/2 ( S ) 

< 1Co(rio)Q{u) + Q{u)+TQ{u) 2 

< (/CoM + iW) 

The last inequality is based on the choice of T as above. The solution r\ is temporally differen- 
tiable to get the equation 

(3.2.12) d t {d t ri) + uidi(d t ri) + u 2 <9 2 (<9tT7) = d t u 3 - d t uidir] - d t u 2 d 2 r] 
So we have the estimate 

( 3 - 2 - 13 )||^i 2r ?|Loo^2iV-5/2( S ) ^ ||m||^ o^2JV-5/2( E ) || d t T] || |oo^2AT-3/2 ( E ) + || O t u\\ ^ R 2N -S/2 ( S ) 

+ \\9tU\\ 2 LXH 2N-5/2^ ||??||^oo H 2iV-3/2( E ) 



< ^/C (r/o) + lJ(^Q(n) 2 + Q(n) 
In a similar fashion, we can achieve the estimate for dfrf as j = 1, . . . , N. 

(3-2.14) dfrj 2 ; < f/C (»yo) + V Q(u) 1 

V ' 11 L oo^2iV-2j+3/2( S ) ~ \^ UV ,Uy y v ' 

To sum up, we have 

(3.2.15) £(r?)<(l + /CoM)P(l + Q(u)) 



Step 3: Estimate of £>(??): 

A straightforward calculation reveals 



(3-2.16) ||r/||| 2if 2Jv+i/2( S ) < T WyWlaetfiN+i/ap) i$ T (^o(vo) + l)j 

Similar to above argument, we have 

(3.2.17) \\d t r]\\ 2 L 2 H 2N-l/2^ < ||n||^ 2/i -2JV-l/2( S ) |M|£cx>#2iV+l/2( S ) + \\u 3 \\ 2 L 2 H 2N-l/2^ 

Furthermore, we have 

112 11^ 2 2 2 

(3.2. 18) 1 1 77 1 1 ^#2^-3/2 ( S ) ^ ||u|| iooH 2iV-3/2( E ) || d t T] \\ L 2 H 2N-l/2 ( S ) + 1 1 <%w| | L 2 H 2N-3/2 ( S ) 

+ 1 1 <9t^ 1 1^2 #2^-3/2(2) 11^11^00^2^-1/2(2) 



< (KqM + lj (^Q(n) 2 + Q(n) 
In the same fashion, we can achieve the estimate for dfrj as j = 1, . . . , N + 1. 

( 3 - 2 - 19 ) lh'HlL-^/2 s fata) + x ) E 2 ^ 

v 7 1=1 
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To sum up, we have 
(3.2.20) 



v( v )<(i + )Co(vo))P(i + Q(u)) 



Step 4: Synthesis: 

It is easy to see the synthesized estimate 



(3.2.21) 



K(r,) = £( V ) + T>(r,) < (1 + /C (%))P(1 + Q(u)) 



□ 



dti + uid\(, + U2d 2 £ = u 3 - uidirjQ - u 2 <9 2 r/o on £ 

e(o) = o 



Next, we introduce a lemma to describe the different between 77 and 770 in a small time period. 

Lemma 3.2. // Q(u) + £0(770) < 00, i/ien /or any e > 0, £/iere exists a T(e) > defending on 
Q(u) and £0(770), swc/i £/ia£ /or any < T < T, we have the estimate 

(3-2.22) ||t7 - 77o||^ ooif 2i V+ i/2 (E) < e 

Proof. £ = 77 — 770 satisfies the transport equation 

(3.2.23) 

By the transport estimate stated above, we have 
(3.2.24) || £|| iooH 2iv+i/2( E ) 

< exp J \\u(t)\\ H 2N+i/2^ dtj ||u 3 (t) - ui(t)dirj Q - u 2 (t)d 2 m\\H 2N + 1 / 2 {S) dt 

< TQ(u)K Q ( m ) 

Hence, when T is small enough, our result naturally follows. □ 

3.3. Forcing Estimates. Now we need to estimate the forcing terms which appears on the 
right-hand side of the linear Navier-Stokes equation. The forcing terms is as follows. 

(3.3.1) F = d t f]bKd 3 u-u-V A u 

(3.3.2) H = r]M 

Recall that we define the forcing quantities as follows. 

N-i „ N 



(3.3.3) K{F,H) = 



L 2 H 2N-2j 



N-l 



(3.3.4) tC (F,H) 



£ j oo^2N-2j-2 



H 2N-2j-2 



3=0 
N-l 



j=0 
2 

L°°_f/2JV-2j-3/2( S ) 



L 2^2iV-2j-l/2( E ) 



E \\ d t F 

3=0 

N-l 

3=0 

N-l 

E H F (°) 

3=0 3=0 

In the estimate of Navier-Stokes-transport system, we also need some other forcing quantities. 

N-l „ N 



E w°) 



H 2AT-2j-3/2( S ) 



(315)J(F,F) = E \ft F 
3=0 

N-l 

(3.3.6) H(F,H) = E|K F 
j"=o 



L 2 H 2N-2j 



L 2 H 2N-2r 



3=0 



L c aH 2N-2j-l/2(Y,) 



N-l 

, + E ||^ 

J"=0 



L2 ff 2JV-23-l/2( S ) 
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Theorem 3.3. The forcing terms satisfy the estimate 
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K(F,H) 



< 



l + K{nf 1 + Q(« 



(3.3.7) 
(3.3.8) 
(3.3.9) 
(3.3.10) 



Proof. The proof is standard, just similar to what we did in derive (2.4.34) to (2.4.47), so we 
omit the details here. In (3.3.7), note the trivial bound 



Hf.h) < ^i + zc^J^i + g^) 2 

H(F,H) < rfl+/Cfa) 2 Vl + Q(u) 2 
)C (F,H) < fl+ZCoM^Vl + ZCoM^ 



(3.3.11) 



L 2 H° 



Also in (3.3.9), the key part is the appearance of bounding constant T. We first trivially bound 
it as follows. 



N-l 

E |$> 

3=0 



L 2 H 2N-2j- 



N-l 
! + E ||^ 

J"=0 

AT-1 

eH^ 

j=0 



lacJj2N-2j 



N-l 

, + E ||^ 

3=0 



£ooJJ2JV-2j-1/2( E ) 



Then just as previous, we can easily show the result. 



□ 



Remark 3.4. The reason why we can get the bounding constant T lies in that, in the nonlinear 
Navier-Stokes equation, the u in nonlinear terms actually has one less derivative than the linear 
part, which makes it available to use lemma A. 14- The appearance of T in (3.3.9) will play a 
key role in the following nonlinear iteration argument. 



4. Navier-Stokes-Transport System 
In this section, we will study the nonlinear system 

dtu — dtfjbKdsu + u • V ' j\u — A_4U + V^p = in Q 

V.4 ■ u = in Q 

Sa(p,u)M = fjjV on S 

u = on Sfc 
u(x, 0) = uq{x) 



(4.1) 



dtV + uidin + «2<9 2 r? = n 3 
, V(x',0) = r/o(x') 



on S 



We will first show a revised version of construction of the initial data and then employ an 
iteration argument to show the wellposedness of this system. 

4.1. Initial Data and Compatible Condition. In the construction of initial data for linear 
Navier-Stokes equation and transport equation, we always assume the other variables are known 
at the moment. However, in the Naiver-Stokes-transport system, (u,p) and n are coupled. That 
means we have to figure out a coupling construction method for this system. 
Define 



(4.1.1) 



/C = /C («o) +K.o(vo) 
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We need a more exact enumeration of terms in T-Lo(u,p), TLoiv) an d JCq(F, H). For j = 0, . . . ,N— 
1, we define 



(4.1.2) 

For j = 1, 
(4.1.3) 

(4.1.4) 



1=0 



JJ2N-21-2 



+ 



d l t H(0) 



H 2N -2I-3/2( 2 ) 



, N, we define 

^■0(V) = ||%||^2JV+l/2(2) 



Ko(v) = IMI^2iv + i/2 ( 2) + \\ d tV(°) 



1=1 



^2JV-2!+3/2(2) 



For j = 1, . . . , N, we write 
(4.1.5) 



(4.1.6) 



Ho{u,p) = \\uo\\%™ 



2 

JJ2N-21 



JJ2N-21-1 



1=0 1=0 

The following lemma records more detailed estimates of above quantities. 

Lemma 4.1. // we write Pj(-, •) /or a polynomial with P(0, 0) = 0. T/ien /or F and H defined 
as in (3.3.1) and j = 0, . . . , N — 1, we /icwe 

(4.1.7) /C J (F,tf) < PjiH^Wrt&P)) 

FordfF(0) andd^Hfi) defined as in (3.3.1) combined with construction in linear problem, and 
j = 0, . . . , N — 1, we have 

(4.1.8) \\Fi(0)\\ 2 H2N _ 2j _ 2 

+ ||^' ? (0)|| j H'2Ar-2j-3/2(2) 

For j = 0, . . . , N, it holds that 



(4.1.9) 
For j = 0, 

(4.1.10) 

Also, 
(4.1.11) 



diu(0)-Diu(0) 



If2N-2j 



<Pj(ni(v),ni(^p)) 



, N — 1, it holds that 
j 



1=0 



^2JV-2j+3/2(2) 

l u - M)||/f2JV-i/2(2) ^5 II^oIIjW (1 + ||?7o||^-2iV+l/2(2)) 



Proof. These bounds are easy to derive by arguing as in lemma 2.31, so we omit the details 
here. □ 

This lemma allows us to combine all the construction of initial data in previous sections under 
the condition that JCq(uq) + /Co(t?o) < °°- Define 9(7/(0) = u$ - Mo, so estimate (4.1.11) implies 
ll^ r ?(0)ll_f/2Jv-i/2(2) ^ ^o- Then we use this for estimate (4.1.7) and have 

(4.1.12) /C°(F,P)<P(/C ) 

for a polynomial P(-) satisfying P(0) = 0. In the following, this polynomial can change from 
line to line. Since we do not care about the detailed form of it, we just omit it. 
Suppose now for j £ [0, N — 2] and that d l t u(0) are known for I = 0, . . . , j, d l t n(0) are known for 
Z = 0,...,j + 1 andd|p(0) are known for I = 0, . . . , j — 1 and that 

(4.1.13) ni(u,p)+ni +1 ( V )+lCl(F,H) < P(/C ) 
According to estimate (4.1.8) and (4.1.9), we know 



(4.1.14) \\Fi(0)\\ 2 H2N _ 2j _ 2 

+ ||^(0)||^2JV-2j-3/2(2) 

+ \\D 3 t u(0) 



< 



ff2N-2j 



P{Ko) 
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By virtue of estimate (2.4.45) to (2.4.47), we have 

2 

H 2N-2j-3 

t>(F>(0),Diu(0)) 



(4.1.15) 



f(Fi(0),Diu(0)) 



+ 



g(W(0),Dlu(0)) 
< P(/Co) 



H 2N~2j-3/2 



H 2N-2j-5/2 



This allows us to define &[p(0) as the solution to (2.3.55) with /, g and h given by f, q and f) 
as above. Then proposition 2.22 with k = 2N implies that 

2 



«0) 



< 



P(JCo) 



(4.1.16) 
Now we can freely define 

(4.1.17) D{ +l u(0) = g°(F j (0),Diu(0),d{p(0)) G H 2N ~ 2j ~ 2 

which satisfies the estimate 



(4.1.18) 

Then we define 



d{ +i u(o) 



< 



H 2N-2(j+l) 



(4.1.19) 

which satisfies the estimate 
(4.1.20) 



1=0 



i + i 
i 



d l t M(0)-di- l+1 u(0) 



H 2N-2(j+2)+S/2 ( S ) ~ v uy 



Combining with above estimates, we have 

(4.1.21) nl + \u,p) + %i +2 (r ] ) +ICl + \F,H) < P(/C ) 
By this iterative method, we can derive that 

(4.1.22) K-\u,p) + <(r?) +IC^- 1 (F,H) < P(/C ) 



Then it remains to define d t ~ p(0) and D^u(0). In order to do this, we need to restate the 
compatible condition for the initial data. We say (uo,f?o) satisfies the N th order compatible 
condition if 



(4.1.23) 



' V^ • (D{u(0)) =0 in n 

D{u(0) = on Ej 

k U o (W(0)+B Ao Diu(0)M o ) = on E 



for j = 0, ...,N-1. 

Then our argument is quite standard, following exactly the path in linear problem, we can 
define df~ 1 p(0) as the weak solution of (2.3.55) with forcing terms defined as in linear problem. 
Then define 



(4.1.24) 



d?u(o) = g^F^io), D?-t u (o),d? -yo)) 



N-l„ 



with corresponding estimates. 

Then we have a final version of initial estimate. 

Theorem 4.2. Suppose («o,t/o) satisfies JCq < oo. Let &[u(0), &[n(0) for j = 0,...,N and 
dfp(0) for j = 0, . . . , N — 1 defined as above. Then 

(4.1.25) /Co < n (u,p) + n (v) < P(Ko) 

for a given polynomial P(-) with -P(O) = 0. 

Proof. Just summarizing all above estimates in the construction will show the result. □ 
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4.2. Construction of Iteration. For given initial data (uo,r]o), by extension lemma A. 15, 
there exists it denned in Q x [0, oo) such that Q(u°) < /Co(^o) achieving the initial data to N th 
order. By solving transport equation with respect to u°, theorem 3.1 implies that there exists 
rj° defined in Q x [0,T ) such that K{rf) < (1 +/C (?fc))(l + Q(u°)) N < K% N . This is our start 
point. 

For any integer m > 1, define the approximate solution (u m ,p m ,n m ) on the existence interval 
[0, T m ) by the following iteration. 

f d t u m - A Am -iu m + V A m-ip m = dtr\ m - x bK m - x d-su m - x in n 

V A m-i ■ u m = in n 

( p ™ 1 _ B^-iu^M 771 - 1 = r f n - 1 M m - 1 on S 

u m = on T, b 



(4.2.1) 



k dtV m + uT diV m + ufd 2 v m = u% on E 

where (u m ,n m ) achieves the same initial data (uo,r/o), and A m , Af m , K m are given in terms of 

This is only a formal definition of iteration. In the following theorems, we will finally prove that 
this approximate sequence can be defined for any m G N and the existence interval T m will not 
shrink to as m — > oo. 

4.3. Boundedness Theorem. 

Theorem 4.3. Assume J° > 5 > and the initial data (uo,?7o) satisfies the N th compatible 
condition. Then there exists a constant < Z < oo and < T < 1 depending on /Co, such that 
ifQ <T <T and ICo < oo, then there exists an infinite sequence (u m ,p m , n m )^ =0 satisfying the 
iteration equation (4-2.1) within the existence interval [0,T) and the following properties. 

(1) The iteration sequence satisfies the estimate 

Q(u m )+JC(n m ) < Z 

for arbitrary m, where the temporal norm is taken with respect to T. 

(2) For any m, J m (t) > 5/2 forO<t<T. 

Remark 4.4. Before we start to prove this boundedness result, it is useful to notice that, 
based on the linear estimate (2.4-54) and forcing estimate (3.3), this sequence can always be 
constructed and we can directly derive an estimate 

(4.3.1) JC(u m+1 ,p m+1 ) + /C(r/ m+1 ) < CP(1 + Q(u m ) + K{n m ) + Kq) exp (tP(1 + JC(n m ))^J 

For a universal constant C > 0. Since the initial data can be arbitrarily large, this estimate 
cannot meet our requirement. Hence, we have to go back to the energy structure and derive a 
stronger estimate. However, this result naturally implies a lemma, which will be used in the 
following: 
If 

(4.3.2) Q{u m - 1 ) + KH 71 - 1 ) < Z 
then 

(4.3.3) IC(u m ,p m )+IC(v m ) <CP(l + ICo + Z) 

Proof. Let's denote the above two assertions related to m as statement P m . We use induction 
to prove this theorem. To note that in the following, we will extensively use the notation P(-), 
however, these polynomials should be understood as explicitly given, but not necessarily written 
out here. Also they can change from line to line. 

Step 1: Po case: 

This is only related to the initial data. Obviously, the construction of u° leads to Q(u°) < /Co- 
By transport estimate (3.1), we have K,{rp) < P(/Cq). We can choose Z > P(/Cq), so the first 
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assertion is verified. 

Define £° = rp — rjo the difference between free surface at later time and its initial data. Then 
the estimate in lemma 3.2 implies ||£°|| LO o#5/2 < TZ. Naturally, sup te [ 0T ] |J°(i) — J°(0)| < 
||C || L oo H 5/2 < TZ. Thus if we take T < 5/(2Z), then J°(t) > 5/2. So the second assertion 
is also verified. In a similar fashion, we can verify the closeness assumption we made in linear 
Navier-Stokes equation, i.e. r\ and ?7o is close enough within [0, T] can also be verified. Hence, 
Po is true. 

In the following, we will assume P m _i is true for m > 1 and prove P m is also true. As long as we 
can show this, by induction, P n is valid for arbitrary n £ N. Certainly, the induction hypothesis 
and above remark implies Q(u m_1 ) + /C« t ~ 1 ) < Z and /C(« m ,p m )+/C(r? m ) < CP(1 + /C + Z). 

Step 2: P m case: estimate of u m via energy structure 

By theorem 2.36, the pair (D^u m ,d^u m ) satisfies the equation (2.4.56) in the weak sense, i.e. 

' d t {D?u m ) - A A m-i (D?u m ) + V^ m -i {d?p m ) = F N in fi 

V^m-i • (D N tu m ) = in n 

S A m-i {D?u m , d?p m )N = H N on S 

Df u m = on S b 



(4.3.4) 



where F N and H N is given in terms of u m and n m 1 . Hence, we have the standard weak 

formulation, i.e. for any test function tp £ X m ~ l , the following holds 

(4.3.5) 

Therefore, when we plug in the test function tp = D^u m , we have the natural energy structure 

(4.3.6) \ ! J \D?u m \ 2 + \f I J \B Am -.D?u m \ 2 = 

i / J(0)|A^ m (0)| 2 + ^ T / a t j|A^ m | 2 + f [ JF N -D?u m - f [ H N -D?u m 
1 Jn 1 Jo Jn Jo Jn Jo Jn 

A preliminary estimate is as follows. 

LHS > \\n?u m \\ 2 LooH0 + \\D?u m \\ 2 L2Hl 
RHS < P(lCo) + TZ\\D?u m \\ 2 LOCHO 

J-\/T"7 II T?N\\ \\r>N m\\ , fifi 1 1 ttN 1 1 1 1 r\N m 1 1 

+ VJ ^lr Hl^ifO || U t U \\ LoaH "I" V 1 \\H 11^00^-1/2(2)11^^ IIl 2 h 1 /2( S ) 

< p(/Co) + || A^lllotfo + Vr ll^ll'^o + v^z 2 1| A^ m ||'oo ff o + \\h n \\ 2 looh _ 1/2{i:) 

+VT\\D?u m \\ 2 L2Hl 

Taking T < 1/(1QZ 4 ) and absorbing the extra term on RHS into LHS implies 

(4.3.7) \\D?u m \\ 2 LaoH0 + \\D?u m \\ 2 L2Hl < P(/C ) + VT\\F n \\ 2 l2h0 + VT ||^||^_ 1/a(E) 

Drop the ||D^ V M m ||^ oo ^ term and we derive further the estimate for d^u m 
(4.3.8) 

^ < P(/C ) +^11^11^0 + ^II^H^-V^E) + II " 

Then we need to estimate each term on RHS. For the middle two terms, it suffices to show it is 
bounded, however, for the last term, we need a temporal constant T within the estimate, which 
can be done by lemma A. 14. Since these estimates is similar to the proof of lemma 2.31, we 
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will not give the details here. 



F 



AMI 2 



H 



iV ||2 



, N-l 
V 3=0 

< P(1+IC + Z)+T(F,H) 

■ N-l 



2 

L 2 H 2 



N-l 

+ E ||^ m 

3=0 



L 2 m 



F(F, H) 



< P(l+ICo + Z) + F(F,H) 

iN mil 



2 

L°°H 2 



J'=0 7 



JV-1 



< TP(l + K(rr- x )) (EP 

< TP(l + /C + £) 

Therefore, to sum up, we have 

(4.3.9) \K um \\l 2 m < P ( /C o) + + /C + -Z) + VTF(F, H) 



2 



(4.3.10) 



Step 3: P m case: estimate of u m via elliptic estimate 

For < n < N — 1, the n t/l order Navier-Stokes equation is as follows. 

' d t {D\ l u m ) - A A (D?u m ) + V_4(di> m ) = F n m n 

• (Df u m ) = in n 

S A {D?u m ,d?p m )N = H n on E 

_ D?u m = on S 6 

where F n and if™ is given in terms of u m and r/™ -1 . 

A straightforward application of elliptic estimate reveals the fact. 

(4.3.11) ||-D™ti m ||^2^2JV-2n+l ^ ||9jZ?™n m ||^2^2JV-2n-l + 1 1 F™ \ | £2 ff2N~2n-l + || H 71 1| jJ2jj2N-2n-l/2 

A more reasonable form is as follows. 
(4.3.12) 

\\an„,m\\2 <■ ||77in||2 _^llH~ n ll 2 j_ II £) n +lo, m II 2 

ll a t U llL 2 /f' 2JV - 2 ™+ 1 ~ lr \\L 2 H 2N - 2n - 1 + II-" llL2^2JV-2n-l/2 + 1 1 C t U || L 2 H 2N-2n-l 

+ ||a t (Z?> m - ^OII^A™-*-! + ||Z?> m - d?U m \\ 2 L 2 H 2N-2n + l 

Then we give a detailed estimate for each term on RHS. These estimates can be easily obtained 
as what we did before, so we omit the details here. It is noticeable that the appearance of T is 
due to lemma A. 14. 



mil^-*-! $ rp(i+/c(r? m - i ))f e k« n 

V i=o 

< TP(l+/C + ^) + ^(F, J ff) 



JV-2 



AT-2 



£oo^2AT-2j 



. + E 



£ j ooJj2N-2j-2 



+ H(F,H) 



I U n ll 2 

l-H || L 2#2iV-2„-l/2( E ) 



V i=o 

< TP(l + /Co + .Z)+H(F,tf) 



Ar-2 



AT-2 



£oo#2JV-2j 



, N-l 

\d t {D?u™-d?u™)\\ 2 L 2 H 2 N -2 n -, < TP(l + /C(r ? — 1 ))( E||^' 



.! + E 



j=0 



£ j oof[2N-2j-2 



+ H(F,H) 



Lootf2N-2j-l 



< TP(l + K, + Z) 
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■ N-2 



f II 2 \ 

V j=o " " / 

< TP(l+K, + Z) 



Therefore, to sum up, we have 

N-l 

(4.3.13) Yl \\d?u m \\ 2 L2H2N - 2n+1 < TP(1 + JC + Z) + H{F, H) 



n=0 

Step 4: P m case: synthesis of estimate for u m 
Combining (4.3.9) and (4.3.13) with lemma A. 12 implies 

(4.3.14) Q{u m ) < P(/C ) + VfP(l + JC + Z) + VTP(P, i?) + ff) 
Then, by forcing estimate (3.3), we have 

HF°,H°) < (l + ^)j(l + Q(«)j<2 4 

H(F°,H°) < r(l+lC( V ) 2 ^(l + Q( U ) 2 ^<TZ 4 
Hence, to sum up, we achieve the estimate 

(4.3.15) Q(u m ) <p()C ) + VTP(l+)Co + Z) 
This is actually 

(4.3.16) Q{u m ) < CP(!C ) + VTCP{\ + K, + Z) 

for some universal constant C > 0. So we can take Z > 2CP{K.q). If T is sufficiently small 
depending on Z, we can bound Q(u m ) < 2CP(K,q) < Z. 

Step 5: P m case: estimate of r] m via transport estimate 
Employing transport estimate in lemma 3.1 

(4.3.17) Kit 1 ) < IC P(Q(u m )) 
Hence, we can take Z = 2JCqP(2CP()Cq)), then we have 

(4.3.18) /C(rT) < Z 



Step 6: P m case: estimate of J m {t) 

Define £ m = r\ m — rjo, then transport estimate in lemma 3.2 implies ||^ m || L00 ^5/2 < TZ. Nat- 
urally, sup tg[0iT] | J m (t) - J m (0)| < ||e m || L oo H 5/ 2 < TZ. Thus if we take T < 5/(2Z), then 
J m (t) > 5/2. A similar argument can justify the closeness assumption of r] m and 770 in studying 
linear Navier-Stokes equation. 

Synthesis: 

Above estimates reveals that if we take Z = CP{K,q) where the polynomial can be given 
explicitly by summarizing all above and T small enough depending on Z, we have 

(4.3.19) Q{u m ) + KH 71 ) < Z 
and 

(4.3.20) J m {t)>5/2 for t € [Q,T\ 

Therefore, P m is verified. By induction, we conclude that P n is valid for any n G N. 

□ 
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Theorem 4.5. Assume exactly the same condition as Theorem (4-3), then actually we have 
the estimate 

(4.3.21) lC(u m ,p m ) + JC{n m ) < P(l + /C + Z) 

Proof. This is a natural corollary of above theorem. Since we know Q{u m ~ l ) < Z, by above 
remark, it implies K,{u m ,p m ) is bounded for any m G N. For convenience, in the following we 
also call this bound Z. □ 

Remark 4.6. /n a&owe proof, the bounding polynomial is changing from line to line, however, 
it can always be explicitly given and does not depend on the data and iterative index m. Since it 
is too complicated and does not help us a lot to understand the problem, we omit its expression 
here. 

AA. Contraction Theorem. Define 

(4 4 1) ^' 9;T ) = IMli-H2 + IMliaffS + \\d t v\\loo H0 + \\dtv\\ 2 L 2 H i + \\q\\ 2 L oo H i + \\q\\ 2 L 2 H 2 

9tt((;T) = ||CHlooH5/2(E) + II^CIIloo H 3/2( S ) + ||^tC||£,2 H l/2( E ) 

Theorem 4.7. For j = 1,2, Suppose that v 3 , q 3 ' , w 3 and Q 3 achieve the same initial condition 
for different j and satisfy 



(4.4.2) 



dtv* - A 



7 a^ j 



V 



A? ■ vJ 







S A ,(q J .' J )X> = C 3 N 3 
v> = 



w 



dtQ 3 bK 3 d^w 3 — w 3 • V A jW 3 in fl 

in Q 

on £ 

on Sft 

in f2 



where A 3 , N 3 and K 3 are in terms of ( 3 . Suppose K.{w 3 ,0), IC(v 3 ,q 3 ) and K.(Q 3 ) is bounded by 
Z. Then there exists < T < 1 such that for any < T < T, we have the following contraction 
relation 



Kiv 1 - v 2 , q 1 - q 2 ;T) < -^(w 1 - w 2 , 0; T) 



(4.4.3) 

(4.4.4) JOT^-C 2 ^)^^ 1 -^^;^ 
Proof. We divide this proof into several steps. 



define v = v 1 — v z , q = q 1 — q l , w = w 1 — w 2 and ( = C 1 — ( 2 , which has trivial initial condition. 



Step 1: Lower order equations 

12 12 

; —v,q = q— q , '. 
Then they satisfy the equation as follows. 

d t v + V^i • S A i (q, v)=H 1 + V A i • I 
V Al -v = H 2 

S Al (q, v)M l =H 3 + D^i.^vW 1 



(4.4.5) 







in $7 
on £ 
on £5 



where 
(4.4.6) 
H 1 = 



H 2 
H 3 



V^i.^QB^v 2 ) - (A 1 - A 2 )Vp 2 + OtCbK^sw + d t CbK l d^w 2 
+d t ( 1 b(K 1 - K 2 )d 3 w 2 - w ■ V^w 1 - w 2 ■ V A iw - w 2 ■ V^i-^to 2 
-V(_4i__42) • V 2 

-q 2 (M 1 - M 2 ) + B A ^v 2 {M l - M 2 ) - B {A i_ A 2 )V 2 {N l - M 2 ) + CAA 1 + C 2 (A^ - M 2 ) 

The solutions are sufficiently regular for us to differentiate in time, which is the following 
equation. 

d t {d t v) + V_4i • S A i (d t q, d t v) = H 1 + V^i • 0®(dtA 1 -d t A 2 ) v2 ) in n 
V Al ■ d t v = H 2 in Q 

S Al {dtq,dtv)J^ =H 3 + B {dtA1 _ dtA 2 ) v 2 Af 1 on S 

k dtv = on 



(4.4.7) 
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where 
(4.4.8) 

H 1 = dtH 1 + Va^i • (B (A i_ A 2 )V 2 ) + V^i • (n {A i_ A 2 } d t v 2 ) + V 9tA i ■ (B A iv) + V A i{B dtA iv) - V 9tA iq 

H 2 = d t H 2 -V dtA i-v 

H 3 = 8 t H 3 + B^.^d^M 1 + U>( A i_ A 2}V 2 d t N 1 - S A i(q, v)8 t M 1 + B^vM 1 



Step 2: Energy evolution for dfV 

Multiply J l dtv on both sides to get the natural energy structure: 

\ [ \dtv\ 2 J l + \ f [ l^d^J 1 = f f J 1 H 1 -d t v+ f f J l H 2 d t q- f f H 3 ■ d t v 
1 Jn 1 Jo Jn Jo Jn Jo Jn Jo Jv. 

+ \ So L idM2dtjl ~Uo L jiB ^ Ai ~ d ^ y2 ■ ° Aidtv 

LHS is simply the energy and dissipation term, so we focus on estimate of RHS. 



^H^dtv < Z 



o Jn 



H 1 



< VtzVw 



H° 

H 1 



dtv\\ H o < Z \\d t v\\ 



L 2 H° 



L 2 H° 



f t J l H 2 d t q < [ J\H 2 - f I (d^qH 2 + j\d t H 2 
Jo Jn Jn Jo Jn 



< 



< 

1-NJ 



/ J 1 qH 2 + \\q\\ LOOH o VfZ 
Jn 



L 2 H° 
+ 



Vtz 



d t H 2 



L 2 H° 



L 2 H° 



H 3 -d t v < 



JE 



j J l qH 2 + VTZVw(^ 

\\dtv\\ Hl/ 2p) < 



dtH 2 



L 2 H0 



H-V2(E) 



H- 



L°°ff- 1 /2(E) 



f I \d t v\ 2 dtJ 1 < ll^H^rii^ill^o 

Jo Jn 

< TZVX" 

So we need several estimate on H l . For the following terms, it suffices to show they are bounded. 
We use the usual way to estimate these quadratic terms. Since we have repeatedly used this 
method, we will not give the details here. 

£ Z ( HCII^S/a + \\d t (\\ L 2 H l/2 + ||^ 2 C|| L2 ^- 1/2 + ||U>|| £ 2 H 1 + \\dtW\\ L 2 H l 
+ IML 2 ^ + |M| L 2 H 1 

< z( V^t" + V^" + Vm' 



H 1 


< A 




L 2 H° ~ \ 



H 2 


< z( 




L 2 H° \ 



d t H 2 



L 2 H° 



< 2^\\d t c\\ L 2 H1/2 + \\a L 2 H1/2 + \\v\\ L 2 H ^ <z(vw+v%t) 

~ Z ( 11^^11^2^1/2 + ||^C|| L 2 H l/2 + ||CllL2 H l/2 + ll^tVll^Hl + \\V\\ L 2 H 1 
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£ooH-l/2(E) 



^ Z { IICIIz,oo ff l/2 + ||5 t C|| L oc^l/2 + \\q\\ L ^ H l + \\v\\ LaoH 2 



< z^Vw + Vt^j 

Also we have the following estimate. 



J l qH 2 < Z\\q\\ LoaH0 



L oo H 



< 



C satisfies the equation 

/ d t c + wjSiC + u^fcC = -N 2 ■ w 
\ c(0) = o 

Employing transport estimate and the boundedness of higher order norms, we have 

||CII LO o H 5 /2 < Vtz\\w\\ l2h5/2 < VTzVw 

dtC satisfies the equation 

/ d t (d t o + wld^dtO + w\d 2 {d t Q = -N 2 ■ d t w - d t N 2 ■ w - a^j^c - d t w\d 2 c 
\ d t c(o) = o 

Similar argument as above shows that 

\\d t c\\ LaoH1/2 < Vtz ( \\d t w\\ L2H1/2 + \\w\\ L2Hl/2 ) < VtzVw 

Combining above transport estimate and lemma A. 12, we have the final version 

/ J l qH 2 < VtzVwVw" + zVw?\\v\\ looh1 

Jn 

< VtzVwVw + zVw( \\v\\ L2H2 + \\d tV \\ L2H o \ 

< VtzVwVw + VfzVw(^ \\v\\ LooH2 + \\d t v\\ LooH0 ) 

Then we consider simplify the last term in RHS of energy structure. 

J J J^idtAi-dtA^ 2 ■ ®AidtV < ^ \\J l ^{d t A^-d t A*)V 2 \\ 2 L2H » + e \pA^tv\\l*H° 

where e should be determined in order for the second term in RHS to be absorbed in LHS. 

P^W-W^vho ^^T\\J ln (d t A^d t A^ 2 \\ 2 L oo H0 < ^TZM 
To summarize, using Cauchy inequality 

(4.4.9) \\d t v\\loo H o + \\d t v\\ 2 L2m < VTz(vi w + m v + m) 



Step 3: Elliptic estimate for v 

based on standard elliptic regularity theory, we have the estimate 

\\d t v\\ 2 Hr + \\H% r + \\H 2 

,2\|| 2 , lln. . . „.2 A/"! II 2 



+ HV41 • (0 ( A^Ai)V 2 )\\ H r + llB^i-^)^ 1 !!^! 



2(E) 
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set r = and take L°° on both sides 

|M|ioo H 2 + IMlioc^l < \\d t v\\ 2 L oo H + H^HiooffO + ||^ 2 ||ioo H l + 1 1 H ^ \ | £oc H \ (S) 

set r = 1 and take L 2 on both sides 

IMl£ 2fl3 + ||g||£ 2Jja < ||^n 2 2 H1 + H^ 1 !!^! + 1|^ 2 ||^2 + ||i? 3 ||^ 2H 3 (E) 

+ HV41 • (0(^-^)^)11^^ + llD^i-^^H^^j 

We need to estimate all the RHS terms. We can employ the usual way to estimate quadratic 
terms in L 2 H k norm, however, for L°°H h norm, we use lemma A. 12. Then we have the estimate 

2 



RHS < \\d t v\\i 2Hl + WdtvWi^ + TZ\W> + m 

To summarize, we have 
(4.4.10) 

IMlLoff' + h\\l-m + \Mhm + WQWhm Z \\dtv\\ 2 L 2 H i + \\dtv\\ 2 LOOH o + tz (V + vnj 

In total of (4.4.9) and (4.4.10), we get the succinct form of estimate 

(4.4.11) w < Vfz(m w + w + m) 



Step 4: Transport estimate for £ 
( satisfies the equation 

/ d t C + w{diC + w\d 2 C = -M 2 ■ w 
\ c(0) = 

A straightforward application of transport estimate can show 

||CIIloo H 5/ 2(e) < exp^C^ ||u-i|| H 3^ ||AA 2 -u;|| H5/2(E) 

< Vtz\\w\\ 2 l2h5/2{j:) < VTzm w 

Similar to the proof of transport estimate in theorem 3.1, we can use the transport equation to 
estimate the higher order derivatives. 

H3tClli°°.ff3/ 2 (E) ~ IK 2 |L°°if5/2(s) W w \\\^mi 2 {T.) + IICIIlooh5/2(i;) \\ w i\\l^h'^/ 2 (y,) ~ 
In the same fashion, we can easily show 

||d 2 C||^ 1/2 < z<sr 

To sum up, we have the estimate 

(4.4.12) wi<m w 



Synthesis: 

In (4.4.11), for T sufficiently small, we can easily absorbed all W term from RHS to LHS and 
replace all Wl with 91™ to achieve 

(4.4.13) W < sfTZW 

Certainly, the smallness of T can guarantee the contraction. 

□ 
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4.5. Local Wellposedness Theorem. Now we can combine theorem 4.5 and 4.7 to produce 
a unique strong solution to the equation system (1.3.11). 

Theorem 4.8. Assume that (u°,n°) satisfies /Co < oo and that the initial data are constructed 
to satisfy the N th compatible condition. Then there exist < To < 1 such that if < T < Tq, 
then there exists a solution triple (u, p, rj) to the Navier-Stokes equation on the time interval 
[0, T] that achieves the initial data and satisfies 

(4.5.1) K(u,p)+K(r))<CP(Ko) 

for a universal constant C > and a polynomial P(-). The solution is unique among functions 
that achieve the initial data and satisfy )C(u,p) + JC(n) < oo. Moreover, rj is such that the 
mapping $ is a C 2N ~ 2 diffeomorphism for each t G [0, T]. 

Proof. We divide the proof into several steps in the following. We first construct an approx- 
imate sequence as in (4.2.1) and assume T is sufficiently small to guarantee the boundedness 
and contraction theorem is valid. 



Step 1: The approximating sequence 

By boundedness theorem, we can construct infinite triple sequence (u m ,p m ,n m ) satisfying the 
boundedness estimate. Hence the uniform bound allow us to take weak and weak* limits, up 
to the extraction of a subsequence. 



(4.5.2) < 



Qj u m 


dfu 


weakly 


in L 2 ([0,T};H 2N - 2 i +1 (n)) for j = 0, 


...,N 


gj u N+l 




weakly 


in X* 




gj u m 


&\u 


weakly 


-*in L°° ( [0, T] ; H 2N ~ 2 i (0) ) for j = 


0,...,N 


cP t p m — 


&> t p 


weakly 


in L 2 ([0,T];H 2N - 2 i{n)) for j = 0, . . 


,N-1 


d\p m — 


dip 


weakly 


- * in L°°([0, T]; iJ 2JV_2j_1 (Jl)) for j 


= 0,...,N 



and 



(4.5.3) { dfT] m -± dfn weakly in L 2 ([0, T]; iJ 2JV_2j+5 / 2 (E)) for j = 0, . . . , N + 1 

weakly-* in L°°([0, T}; H 2N+l / 2 {^)) 
d{n weakly-* in L°°([0, T]; iJ 2iV_2j+3 / 2 (S)) for j = 0, . . . , N 

According to the weak and weak* lower semi-continuity of the norm, we have that for the limit 
function 





V 


weakly 


Q tV m 


->> d t n 


weakly 


d{n m 


dfn 


weakly 


n ™ 


V 


weakly 


&\n m 


-din 


weakly 



(4.5.4) 



K(u,p)+K( V )<P(]Co) 



Furthermore, since the collection of triple (v, q, Q is closed under the norms rio(u,p) and Ho(n) 
for initial data and (u m ,p m ,n m ) always achieves the same initial data, we deduce that the limit 
(n, p, n) also achieves this. 



Step 2: Contraction 

Now we want to improve the result from weak convergence to strong convergence in the norm 

_ „.m+l 



^JyR{n : T) + 0T(n,p, T). In the contraction theorem, for m > 1, set v 1 



,m+2 



W 1 = u m+1 , W 2 



P 



,m+2 



P 



m+1 



7] 



m+1 



_ Q ur cons t ruc tion satisfies 



all the requirements of contraction theorem, so we can naturally deduce that 



(4.5.5) 

and 
(4.5.6) 



m+2 



u ™+l. p m+2 



P m+1 ;T) < =-m(u 



1 



m+1 



u m. p m+l 



p m ;T) 



Wl(n 



m+1 



n m ;T)<m(u 



m+1 



u m ;p m+1 



p m -T) 
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These estimates implies (u m ,p m ) is Cauchy under the norm \/0t(-, •; T) and n m is Cauchy under 
the norm y/%R(-\T). So as m — > oo, we have 

' u m ^u in L 2 ([0,T];H 3 (n)r\L°°([0,T};H 2 (n)) 

d t u m -»• <%u in L 2 ([0,T];ii' 1 (fi) n L°°([0,T]; 

p m in L 2 ([0,T];H 2 {tl) n L°° ([0,T]; H 1 (0,)) 

rT^n mL°°([0,T];F 5 / 2 (S)) 

d t r) m ^d t r) in L°°([0, T]; # 3 / 2 (£)) 

[ d 2 t n m ^d 2 t n inL 2 ( [0, T] ; */ 2 (£)) 



(4.5.7) 



Step 3: Interpolation 

Since the triple (u m ,p m , rj m ) is bounded under higher order norms and strongly convergent under 
lower order norms, we can use an interpolation argument to get better convergence properties 
of it. The standard interpolation theory implies the following relation. 



(4.5.8) 



dtf 



L 2 H° 



< C(T) 



L 2 H° 



dlf 



L 2 H° 



for j > k > 0, 6 = 6(j, k) G (0, 1) and C(T) a constant depending on T, which reveals that 



(4.5.9) 



d{u m 
df-p m 
cP t n m 



cP t u in L 2 ([0,T];H°(n) for j = 0, . . . , TV - 1 
^p in L 2 ([0,T};H°(n) for j = 0, . . . , N - 2 
d\r\ in L 2 ([0, T];H°(E) for j = 0,...,N 



Then we use spacial interpolation between H° and H k to deduce that 



(4.5.10) 



dip 
V 

dtV 



in 



&iu in L 2 ([0,T];H 2N - 2 i(Q) for j 



0,...,N-l 



■rn 



L 2 ([0, T]]H 2N ~ 2 i~ l (fi) for j = 0, . . . , N - 2 



->• c^?p in 

r? in L 2 ([0,T];H 2N (Z) 

m ^8 t n in L 2 ([0,T];i? 2Ar - 1 /2(S) 

m 



&\n m ^8P t n in L 2 ([0, T]; tf 2A ^'+ 2 (£) for j 
By lemma A. 12, we have the convergence in L°° and C° norm. 



(4.5.11) 



Qj p n 

r/ m - 
d t rf 
4rf 



-)• in ^([O,^;^- 2 ^- 1 ^) /or j 
-> $p in C7 ([0,T] ;j ff 2JV -^- 2 (O) for j 
n mC°([0,7l;fr 2JV - 1 /2(E) 
in C°([0,T];# 2JV - 3 / 2 (£) 
-»• ^'n in C°([0, T]; iJ 2JV_2j+1 (S) /or j 



0, 



0, . 
0,. 



, N 



N -2 
JV-3 



0,...,AT-1 



Step 4: Passing to the limit 

The strong convergence result as above are more than sufficient to pass to the limit in the equa- 
tion (4.2.1). Naturally, the limit triple (u,p,rj) is the strong solution to the problem (1.3.11). 

Step 5: Uniqueness 

Suppose there are two solutions {u l ,p l ,rf) for i = 1,2, then by contraction theorem, we have 

(4.5.12) miu 1 - u 2 y -p 2 ;T) < ^(u 1 - u 2 ,p x - p 2 ;T) 
and 

(4.5.13) 9JT(r/ - r, 1 ; T) < ^(u 1 - u 2 ,p l - p 2 ;T) 
which implies u 1 = u 2 , p 1 = p 2 and n 1 = rj 2 . 
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Step 6: Diffeomorphism 

Since J(t) > 5/2 > 0, so it is sufficient to guarantee diffeomorphism. The smoothness of r\ 
imphes it is in fact a C 2N ~ 2 diffeomorphism. 

□ 



Appendix A. Analytic Tools 

A.l. Products in Sobolev Space. We will need some estimates of the products of functions 
in Sobolev spaces. Since these results have been proved in lemma A.l and lemma A. 2 of [1], we 
will present the statement of the lemmas here without proof. 

Lemma A.l. Let U denote either £ or fi. 

(1) Let < r < 8i < s 2 be such that Sl > n/2. Let f G H Sl (U), g G H S2 (U). Then 
fg G H r {U) and 

(A.l.l) \\fg\\ H r < \\f\\ Hsl \\g\\ H . a 

(2) Let < r < Sl < s 2 be such that s 2 > r + n/2. Let f G H Sl (U), g G H S2 (U). Then 
fg G H r (U) and 

(A-1.2) \\fg\\ Hr < \\f\\ HS1 \\g\\ HS2 

(3) Let < r < Sl < s 2 be such that s 2 > r + n/2. Let f G ^ r (S), g G # S2 (£). Then 
fg G and 

(A.1.3) II/^IIh-i < \\f\\ H -r \\9\\h-2 

Lemma A.2. Suppose that f G C 1 (S) and g G H 1 / 2 ^). Then fg G H 1 / 2 ^) and 
(A.1.4) 11/511^/2 < ll/llci II^IIhi/2 

A. 2. Identities involving ^4. We now record some identities involving A. Since they are 
direct consequences of the definition of A, we won't give detailed proof here. Lemma A. 3 in [1] 
shows the similar results. 

Lemma A. 3. The following hold. 

(1) For each j = 1, 2, 3, we have that dk(JAjk) = 0. 

(2) On S we have that JAe% = N , while on S;, we have that JAe% = v. 

(3) Let R defined by (24.9). Then R T M = -d t M on S. 

A. 3. Poisson Integral in Infinite Case. For a function / defined on S = R 2 , the parame- 
terized Poisson integral is defined by 

(A.3.1) V e f(x',x 3 )= f /(Oe e Kl a *e w *d£ 

where /(£) denotes the Fourier transform of f{x') in i? 2 and < e < 1 is the parameter. 
Although V e f is defined on R 2 x (— oo,0), we will only consider the part R 2 x (— &o>0) here. 

Lemma A. 4. Let V e f be the parameterized Poisson integral of function f which is in homoge- 
neous Sobolev space H q ~ 1 ^ 2 (T,) for ggN. Then we have 

(A.3.2) ||V^/ll^o < j \\f\\% q - 1/2 



where C > is a constant independent of e. Ln particular, we have 
(A.3.3) H^/ll^ < 7 II/II^-i/ 2( e) 
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Proof. By definition of Fourier transform, Fubini theorem and Parseval identity, we may bound 



||V«7>7IIho < {^) 2q / / \Z\ 2q HO e 2 ^ x 

JR 2 J -bit 

n f o - 2 (\ - e~ 2eb °^ 
= (2vr) 2 « / |£| 2 « f(0 ' 

Jr 2 



M£dx 3 < (2vr) 2 « f |e| 29 HO V f e 2 ^dx 
Jr 2 J-bo 



(27T) 



2q 



2e Jr 2 



iei 29_1 HO ^ 



7T 

< - 



, - .Iif9-V2(E) 

We can simply take C = 7r to achieve the estimate. Furthermore, considering 



1 _ e -26& |£| 



< min 



1 



2eier 



and < e < 1, we have 

(A.3.4) \\V*f\\ 2 H0 < ^ H/I&op) 

Hence, (A. 3. 3) easily follows. 

We also need the L°° estimate for Poisson integral. 



□ 



Lemma A. 5. Let V e f be the parameterized Poisson integral of function f. Then for q G N and 
s > 1, we have 



(A.3.5) 



|V^7llLo<CJ/||^ +s(s) 



Proof. A simple application of Sobolev embedding and lemma A. 4 reveals 
||V 9 7>7ll£~ < C\\V^f\\ 2 Hs+1/2 < c ii/n 



□ 



The following lemma illustrate the specialty of derivative in vertical direction. 
Lemma A. 6. Let V e f be the parameterized Poisson integral of function f . Then we have 
(A.3.6) \\d3^f\\ 2 Loa <Ce\\f\\ 2 H5/2(s) 
where C > is a constant independent of e. 
Proof. We can simply bound as follows. 

(A.3.7) llftPVIll- < C \\d 3 V*f\\ 2 H2 < Ce 2 ||P7H 2 H3 < Ce ||/fe 6/a(E) 

Note that the first inequality is based on Sobolev embedding theorem and the third one is a 
straightforward application of lemma A.4, so we only need to verify the second inequality in 
detail. 

By definition, it is easy to see 



(A.3.8) 

Hence, we have 
llftPVI&o = 

< 



d 3 V e f 



= e[ \Z\HO 
Jr 2 



I L ^ M =M (2 * )2 l L if|2 \ m 



-b 

(4 Pir ' ff "" 



d£dx 3 



Similarly, we can show that for i,j = 1, 2, 3 

ll^7ll^o < 



1^711^0 < 



(2vr) 4 

e 2 
(2vr) 6 



\duV*f\\ 2 H0 



\diur £ f\\ 2 H0 
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\\d 3 f\\m<Ce 2 \m 2 H3 



□ 



A. 4. Poisson Integral in Periodic Case. Suppose £ = (L\T) x (L2T), we define the param- 
eterized Poisson integral as follows. 

(A.4.1) r £ f(x',x 3 )= Yl e 27vin - x 'e e ^ X3 f(n) 

n&{L^Z)x{L^ 1 Z) 

where 



(A.4.2) 



f(n) = / f( x >)——dx> 



Lemma A. 7. Let V e f be the parameterized Poisson integral of function f which is in homoge- 
neous Sobolev space 

#?-i/2( S ) f or q G N _ Then we have 



(A.4.3) 



Iv^/IIho < 



c 



J/9-1/2 



where C > is a constant independent of e. In particular, we have 
(AAA) 



\V*f\\ 2 Hq < £ ||/||L- 1/2(S) 



Proof. By Fubini theorem and Parseval identity, we may bound 



|V^/ll^o < (2vr) 2 « £ 

ne(L^ 1 Z)x(L^ 1 2 

= w 2q E 



< 



(27T 



|2g 



2e 



E 



r?. 



n 





-60 

2q 



n\ 2q e 2e ^ X3 



f(n) 



dx 3 



fin) 



2 f\ - e -2efe |n| 



2e|n| 



25-1 



f(n) 



< - 



* II x\\2 



Hi- 1 / 2 



ne(L- 1 Z)x(L- 1 Z) 

We can simply take C = tt to achieve the estimate. Furthermore, considering 



I _ g-2e&oM 

2e\n\ 



< min 



1 



2e|n 



■,6c 



W'fWlo < ^ ||/|||o (s) 



and < e < 1, we have 
(A.4.5) 

Hence, (A.4.4) easily follows. 

Similar to infinite case, we give an L°° estimate for Poisson integral. 



□ 



Lemma A. 8. Let V e f be the parameterized Poisson integral of function f . Then for gGN and 
s > 1, we have 

iiv«p e /ni<»<cii/ii^ + . (E) 

Proof. A simple application of Sobolev embedding and lemma (A. 7) reveals 



(A.4.6) 



|v*7>7lll» < c\\v^f\\ 2 HS+1/2 < c 11/11 



□ 



We still need the following lemma to illustrate the specialty of derivative in vertical direction. 



LOCAL WELL-POSEDNESS OF VISCOUS SURFACE WAVE WITHOUT SURFACE TENSION 47 

Lemma A. 9. Let V e f be the parameterized Poisson integral of function f . Then we have 
(A.4.7) \\d 3 P e f\\ 2 L „<Ce H/ll 

where C > is a constant independent of e. 
Proof. We can simply bound as follows. 

(A.4.8) \\d 3 P e f\\ 2 Lao < C ||<9 3 P7llff2 < Ce 2 \\P*ff H3 < Ce ||/||^/ 2(s) 

Note that the first inequality is based on Sobolev embedding theorem and the third one is a 
straightforward application of lemma (A. 4), so we only need to verify the second inequality in 
details. 

By definition, it is easy to see 

(A.4.9) d 3 V e f = e J2 \n\ f(0e e ^ X3 e 27rix '- n 

ne(L- 1 Z)x(L- 1 Z) 

Hence, we have 



\\dsV*ff H0 = e 2 £ |° \n\ 2 \f(n) 



2 e 26|n|*3 dx3 



(2vr) 2 

J2 



((2vr) 2 f |n| 2 |/(n)| 2 e 2£ H^ dx; 

ne(L- 1 Z)x(L- 1 Z) 



< 77^2 W d ^f Who 



(2vr) 2 

Similarly, we can show that for i, j = 1, 2, 3 

11^71^0 < \\dnV e f\\ 2 H0 

11^711^0 < J^y- ||W7II^0 

So we have 

(A.4.10) \\d3V e \\ 2 H 2 <Ce 2 \\fj% 3 



□ 



A. 5. Poincare-Type Inequality. We need several Poincare-type inequality in fl. Since all 
these lemmas have be proved in lemma A. 13- A. 14 in [1], we will only give the statement here 
without proof. 

Lemma A. 10. It holds that \\f\\ H0 < \\f\\ H i < \\Vf\\ H o for all f £ H 1 ^) such that f = on 

We also need one type of Korn's inequality for both infinite case and periodic case. 
Lemma A. 11. If u = on then \\u\\ H i < ||Du||^o- 

A. 6. Continuity and Temporal Derivative. In the following, we give two important lemmas 
to connect L 2 H k norm and L°°H k norm. 

Lemma A.12. suppose thatue L 2 ([0, T];H S1 (ft)) and d t u G L 2 ([0, T}; H S2 (tt)) for si > s 2 > 
and s = (si + S2)/2. Then u G C°([0, T]; H S (Q)) and satisfies the estimate 

(A.6.1) \\u\\l~hs < II«(0)IIh- + IMlWi + \\dM\ 2 L 2 H s 2 

where the L 2 H k norm and L°°H k norm are evaluated in [0, T] . 
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Proof. Considering the extension theorem in Sobolev space, we only need to prove this result 
in the R n case. The periodic case can be derived in a similar fashion. Using Fourier transform, 

d t \\u(t)f Hs = 2m( [ (0 2s u(tt)dM^)dA <2 / {0 2s m,t)\\dtu^tM 

V JR" J JR n 

JR n JR n JR n 

= IN*)llH'i+IIMt)llH.2 

So integrate with respect to time on [0, t] 

\\u(t)\\ 2 H s < \\u(0)\\ 2 Hs + \\u\\ 2 L2Hsl + \\d t u\\ 2 L2Ha2 

□ 

In [1], lemma A. 4 provides us another form of estimate, which will be employed to prove the 
relation between the topologies generated by these norms. 

Lemma A. 13. Under exactly the same condition as above lemma, we have the estimate 

1 



(A.6.2) 



\L ca H s 



< 1 + 



U 



+ \\dtu\\i 2HS 



The following lemma shows the estimate in another direction. 

Lemma A. 14. For any u G L°°([0,T];H k ) within [0,T], we must have u G L 2 ([0, T];H k ) and 
satisfies the estimate 



(A.6.3) Ihll^^fc < TUnll^^fc 

Proof. The result is simply based on the definition of these two norms. 



□ 



A. 7. Extension Theorem. The following are two extension theorems which will be used to 
construct start point of iteration from initial data in proving wellposedness of Naiver-Stokes- 
transport system. 

Lemma A. 15. Suppose that dju(0) G H 2N ~ 2 i(Vi) for j = 0, . . . , N, then there exists a extension 
u achieving the initial data, such that 

d{u G L 2 ([0,oo);F 2iV - 2j ' +1 (^))nL oo ([0,oo);F 2iV " 2j '(^)) 

for j = 0, . . . , N. Moreover, 



N 



(A.7.1) 
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L 2 H 2N-2j + l 
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&> t u 



L oo H 2N-2j 
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<E||^(°) 

j=0 



H 2N-2j 



Proof. The same as lemma A. 5 in [1]. 



□ 



Lemma A. 16. Suppose that d[p(0) G H 2N 2jf for j = 0, . . . , N — 1, then there exists a 

extension p achieving the initial data, such that 

&>p G L 2 ([0, 00); H 2N ~ 2 :> (n)) n L°°([0, 00); H 2N - 2 i-\tt)) 

for j = 0, . . . , N — 1. Moreover, 



(A.7.2) 



N-l 
3=0 



+ 



L 2 H 2N-2i 

Proof. The same as lemma A. 6 in [1] 
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L oo H 2N-2j 



N-l 
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□ 
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Appendix B. Estimates for Fundamental Equations 
B.l. Transport Estimates. Let £ be either infinite or periodic. Consider the equation 

/ dtV + u ■ Dr) = g in S x (0, T) 

{ ' \ v(t = o) = m 

We have the following estimate of the regularity solution to this equation, which is a particular 
case of a more general result proved in proposition 2.1 of [3]. Note that the result in [3] is stated 
for £ = R 2 , but the same result holds in the periodic case as described in [7]. 

Lemma B.l. Let rj be a solution to equation (B.l.l). Then there exists a universal constant 
C > such that for any s > 3 

(B.l. 2) IMIlooh. <exp fcjf \\Du(r)\\ H 3 / 2 6x\(\\rio\\ H ,+ J \\g(r)\\ H , 6ij 

Proof. Use p = P2 = 2, iV = 2 and £ = s in proposition 2.1 of [7], along with the embedding 

ff3/2 ^ B 2,oo nLO °' D 

B.2. Elliptic Estimates. We need the following elliptic estimate for Stokes equation. 
Lemma B.2. Suppose (u,p) solves the equation 

-Au + Vp = 4> G H r - 2 (Q) in n 



(B.2.1) 

Then for r > 2, 



V • u = ?/> G H r ~ 1 (ft) in n 

(pi - Ou)e 3 = <p € F r - 3 / 2 (£) on £ 
u = on £5 



\\ U \\ht + lblli^-1 < ||0||ffr-2 + llV'llffr-l + ||^||^r-3/2 

Proof. The same as lemma A. 15 in [1]. □ 
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